SHADOWING MATCHING ERRORS FOR 
WAVE-FRONT-LIKE SOLUTIONS 
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Abstract. Consider a singularly perturbed system 

eut = e'^Uxx + f{u, x,e), u e M." , x e R,t > 0. 

Assume that the system has a sequence of regular and internal 
layers occurring alternatively along the a:-direction. These "mul- 
tiple wave" solutions can formally be constructed by matched as- 
ymptotic expansions. To obtain a genuine solution, we derive a 
Spatial Shadowing Lemma which assures the existence of an ex- 
act solution that is near the formal asymptotic series provided 
(1) the residual errors are small in all the layers, and (2) the 
matching errors are small along the lateral boundaries of the adja- 
cent layers. The method should work on some other systems like 



1. Introduction 

One of the most effective method to study wave-front-hke solutions 
in singularly perturbed reaction-diffusion equations is the formal as- 
ymptotic method. Fife [13] studied the formation of sharp wave fronts 
and obtained conditions on the stability of such solutions by the formal 
method. Matched expansions for systems of any number of equations to 
any order in e has also been constructed, [25]. However, the asymptotic 
method does not guarantee that there is an exact solution near the for- 
mal solution. Many efforts have been made to derive methods that ver- 
ify the vahdity of formal solutions, [1, 2, 3, 5, 16, 17, 18, 20, 22, 23, 24]. 
We present a new method in this paper. 

Although the method should work for some higher order parabolic 
systems, for simplicity, consider a second order singularly perturbed 
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system, 

eut^e^u^x + f{u,x,e), e R", x e R, t > 0. (1.1) 

Assume that (1.1) has an infinite sequence of regular and internal 
layers occurring alternatively along the x direction. Assuming by a 
matched expansion, we have constructed the following formal series for 

m 

(1) positions of moving wave fronts r]^{t, e) = X] ^^Vjit): ^ G ^; 



m 

(2) solutions in the £ih regular layer u^{x, t, e) = X! ^''uf\x, t); 



(3) solutions in the £th internal layer u^^{^, e) = ^ e^Uj\^, t), 



where ^ — {x — ri^{t, e))/e. 

Here the superscript "R" and "S" stand for regular and singular (in- 
ternal) layers. 

Although our method can be applied to systems with various kinds 
of boundary conditions. To simplify the illustration, assume that the 
nonlinear term is periodic in x, and we are looking for solutions that 
are periodic in x. Assume that the formal expansions are periodic in 
£, compatible with the period of /. That is, there exists a constant 
Xp > and an integer ip > 0, such that the following holds: 

Periodicity Hypotheses 

1. f{u,X + Xp,€) = f{u,x,€); 

2. r/^+^p(t,e) = r]^{t,e)+Xp; 

3. M-^(^+^p)(x,t,e) =u^^ix -Xp,t,e)] 

4. M^(^+^-)(e,t,e) = M^^(e,t,e). 

Due to the periodicity of / and the formal solutions, all the estimates 
in this paper are uniformly valid with respect to the layer index £, and 
this will not be repeated in the paper. The periodicity plays no other 
rule apart from mentioned above. 

Through out this paper, let < p < 1 be a fixed constant. Let the 
width of the internal layer be 2e*^ Let 

A family of curves F* = {{x,t) : x = y\t),t > 0}, i E Z, divides the 
domain {x G M} x {t > 0} into infinitely many strips, = {{x,t) : 
y^~^{t) < X < y'^{t),t > 0}, i e Z. is an internal (or regular) layer if 
i = 2£ (or i = 2£-l). 
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A doubly infinite sequence will be denoted by {/i*}, with the norm 

II^^IU* = supjll/i^lgi} if each is in a Banach space A sequence of 
functions {w'^{x, t)}, each defined in E*, is called a formal approximation 
subordinate to the partition, if the residual error in E', 

-g' = ewl-e^wl^-f{w\x,e), (1.3) 

and the jump error at the boundary r\ 

i+l 




w^+'), (1.4) 

are small. We also call the piecewise smooth function w which is equal 
to in E* a formal approximation. In singular perturbation prob- 
lems, the formal approximation can naturally be provided by matched 
expansions. For (x, t) e E', 

, i«^(x,t,e), i = 2£-l. 



We also assume that y\t) and w^{t) are defined for all i > and 
approaches y\oo) and w\oo) as i — > oo. (Our method can show the 
existence of exact solutions in finite time if and exist in finite 
time.) 

Under certain conditions, we want to show that the smallness of 
residual and jump errors implies that for each initial condition u{x, 0, e) 
near w{x, 0, e), there is a unique exact solution u near w. To this end, 
let M = -u;* + -u* in E*. The correction terms satisfy nonlinear 

equations with the forcing terms g\ as in (1.3), in E', 

eui = + f{w' + u\ X, e) - f{w\ x, e) + g\ (1.6) 

and jump conditions 5*, as in (1.4), at P, 



{u' - u'+^) = 5\ (1.7) 



The small solutions {u^} can be solved by linearization and contrac- 
tion mappings in certain Banach spaces. The process is much like a 
shadowing lemma argument used in [23]. However, since the jumps 
do not occur in the temporal direction, we can not directly use the 
normal Temporal Shadowing Lemma as in that paper. A new Spatial 
Shadowing Lemma is needed to correct jump errors along the spatial 
direction. 

Several coordinate systems are used in this paper. We denote the 
original coordinate system, as used in (1.1), by 9^(0). Other coordinates 
will be defined through ^(0) and will be denoted by $H(1), ^(2), etc. 
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We will denote a two dimensional region fl {R x /}, where / is a 
time interval, by fl / for simplicity. 

We now illustrates basic ideas used in this paper. Consider a coor- 
dinate change in the region fl [0, eAr], 

fH(0) ^ : n [0, eAr] 9 {x, t) r) G x I^, 

y = x- E{x,t,e), 

C-[y-{y'-' + y')]/e, r^t/e. 

Here y = x — t,e) is a near identity change of coordinates that 
straightens F*"^ & F* so that the functions y^~^ & y^ are time inde- 
pendent. = {-L'{e) < i < L\e)}, L\e) = {y' - y'-^)/2e and 
7, = [0, Ar]. 

Using the coordinate system in the region x we linearize 

system (1.6) at a fixed time r = 0. We are led to 

< = 4^ + y^(04 + ^^(0^^ + e, T, e) + g\ (1.8) 

where is a scalar multipher, and jV* is a small term. Eqn (1.8) is to 
be solved with an initial condition 

<(O=^(e,0,e)-«;^(e,0,e). (1.9) 

The change of variable can also be made such that it does not affect 
the jump condition (1.7) at F*. For illustration, let jV' = 0. Assuming 

that the coefficients of Eqn. (1.8) can be extended to ^ G M, as well 
as Uq and Using the variation of constant formula in a suitable 
Banach space, (1.8) can be solved for ^ G M with nonzero Uq and 
but no boundary conditions. We then only have to treat (1.8) with 

= 0, (7* = and some jump boundary conditions. 

Define H^-^^'^'^-^^Q = H^'^'^ilr) x H^-^\I^). Let u' G H^^\Q} x 

(uH^ .)\ 

Ir), from the Trace Theorem [28], the mapping ^ — > ( ^ 1; ~^ 

^o.75xo.25^j^^ is continuous. However, Eqn (1.8) does not generate 
a flow in H^-'^^''°-^^{Ir). Following the idea of [21, 30], we want to 
find subspaces W%^) © 1^"(0 = H^-^^''^-'^^{Ir) such that the linear 
homogeneous equation associate to (1.8) can be solved for ^ > if 

the boundary value {^^^^ = (^jj G W^{i,) (or ^ < 6 if (^j) G 

VF"(^o))- We also want to show that these solutions decay as ^ moves 
away from ,^o- Then the idea of Temporal Shadowing Lemma can be 
used to solve systems (1.8) and (1.7). 

The dichotomy splitting can easily be achieved by applying the Fourier- 
Laplace transform to (1.8). Using A for images of Laplace transforms. 
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we have a first order system, which will be called a dual system asso- 
ciated to (1.8). 

/ I \ (u\s)\ , / \ 

- [sI-A^iO -V'iOl) [ns)) ^ 1-^' + ^7(1.10) 

Equation (1-10) is an ODE in a Banach space of homomorphic func- 
tions (Hardy-Lebesgue classes). It can be shown that in certain region 
of s, (1.10) has an exponential dichotomy. Then the inverse Laplace 
transform gives the desired dichotomy in i?o'^^^°"^^(/^). 

We have treated the problem in a short time interval. In the coordi- 
nates 9^(0), let Ij = \jAt, {j + l)At], < j < r - 1, and Joo = [tf, oo). 
Here At = eAr and tf = rAt. The constant tf has to be so large 
and At > so small that the variations of y''{t) and w\t) with respect 
to t are small in each time interval. We can then use a near identity 
change of coordinates to straighten the boundary P, i £ Z, and use a 
linearization at a frozen time with very small errors. Nonlinear systems 
like (1.8) can be solved recursively in time intervals Ij, < j < r — 1, 
and 7oo) using the values of at the end of the Ij as the initial data 
for Ij+i- We remark that the number of intervals r-|-l— >ooase— >0. 
The main difficulty in this paper is to control the growth of as r gets 
larger. 

The outline of this paper is as follows. We state our main result, 
the Spatial Shadowing Lemma in §2. The proof of that lemma is also 
included there, but should be read after §6, since Theorems 5.1, 6.1 
and 6.2 are used in the proof. The theory of dichotomy splittings is 
developed in §3 through several technical lemmas. A useful property 
on Evans function, due to Gardner and Jones [19], is used to prove 
Lemma 3.10. Readers who want to know the main flow of proofs can 
skip §3 on the first reading, but come back for those lemmas when they 
are quoted. In §4, we introduce the change of coordinates that is used 
to straighten the curve T\ We also study a system of linear equations 
with jump conditions along their common boundaries (Lemma 4.1). 
The result there is the main tool to handle jump boundary conditions in 
the next two sections. In §5, we treat equations in the final interval Ir^. 
The result obtained there also tells us how small the upper bound of 
Wexact — w\ must be at time t = tf m. order to construct exact solutions 
in /(X). In §6, we treat equations on each time interval Ij, < j < 
r — 1. The error on each Ij has to be small so that the accumulation 
on all r intervals yields a small error at t = tf. In §7, the Spatial 
Shadowing Lemma is applied to the singularly perturbed system (1.1). 
We briefly review the construction of matched asymptotic expansions 
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and formal approximations for (1.1) and the hypotheses used for such 
constructions, [25]. Wc show that these formal approximations satisfy 
the hypotheses of the Spatial Shadowing Lemma, as stated in §2. A 
precise relation between eigenvalues and wave speeds of internal layers 
is given in Lemma 7.1, which may be of some independent interest. 

We have concentrated on global wave-front-hke solutions in this pa- 
per. However, our method can treat short time solutions with much 
weaker hypotheses. These local solutions are discussed in [27]. The re- 
sult obtained there applies to various general systems including Cahn- 
Hilliard eqautions and viscous profile for conservation laws, which do 
not satify our hypotheses in §2. 

2. Main Result 

Let Q and / be space and time intervals respectively. Define the 
following Banach spaces and norms. 

H^'^{n X I) = {u{x,t) -.Qx I \ u, u^^ and Ut G L^{n x 7;R")}. 

\u\H^,i{nxI) = |li|L2 + \Uxx\l2 + 

H''{I) = W'^il; W),r> 0, is the usual Soblev space. 
^rx.Jj^ = X H'{I), r > 0, s > 0. 

For a constant 7 e R, let 

L2(fi X M+,7) = {m : 1] X M+ ^ I e^^^u E L^{VL x M+)}. 

|li|L2(s7xM+,7) = |e~'^*'u|l,2(f2xR+)- 

X R+, 7) = {k : Q X R+ ^ I e-^*u e x R+)}. 

k|j/2,i(f2xM+,7) = |e~'^*Xi|if2,i(nxR+)- 

For a constant 7 < 0, let 

X{VL X R+, 7) = {m : M = Ml + M2, ui G L'^{VL), U2 e L'^{fl x R+, 7)}. 

|'«^U(7) = |mi|l2(q) + |-U2|l2(QxR+,7)- 

X'^'\Q X R+,7) = {u:u = ui + U2,uie H\n), U2 G H^'\n x R+,7)} 

^1x2,1(7) = \ui\h2(Q) + |xi2|if2,i(nxR+,7)- 

It can be verified that ior u e X{Q xR+ , 7) or X'^^^{Q x R+, 7), the 
decomposition u = U1+U2 is unique. We often use simplified notations, 
for example X^'^(7), instead of X'^'^{fl x R^,7). 

Let S = {{x,t) : y\t) <x < y'^{t), t > 0}, where y^, y^ are smooth 
functions of t, be a two dimensional region unbonded in the t direction. 
We say that u e L2(E,7), H^'\i:,-f), X(E,7) or X^'\^,-f), etc., if m 
is the restriction of a function U e L^(R x R+, 7), etc., to the domain 
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E. The norms are defined by 

I'"|l2{E,7) = inf{|?7|L2(RxR+,7)}, 
|m|x2.1(S,7) = illf{|^U2.i(RxR+,7)}- 

Let 

X*^(7) ^{u\u^ui + U2,uie W, U2 e H'^i^)}, 7 < 0. 
Xkixk^^^^ = X''^{-f) X > 0, A;2 > 0. 

If / is a finite time interval, similar definitions can be given to spaces 
L2(fi X /,7), H^'\n X /,7), ii-^.i^j] n j^^)^ but not X(0 x 1,7) or 

As in §1, let {tf *} be a formal approximation for (1.1). Let {rj^} be 
a formal approximation for the wave fronts. We extend the domain of 
w' from EUo R X R+ by letting 



w'^{x, t, e) = 



w\y'-\t),t,e), ifx<y'-\t), 
y{f{t),t,e), iix>f{t). 
Assume 

H 1. There exist C, 7 > such that for all small e and i, i E 1^, 

\w\x,t,e) -w''{x,oo,e)\ < Ce^^*, {x,t) G 
\r]^{t, e) - r/^(oo, e) | + \Dtr]^{t, e) | < Ce'^', t e M+. 
Here w*(a;, 00, e) = lim^^oo w\x, t, e), etc. 

H 2. There exists a > such that in each regular layer E*, i = 2£ — 1, 

Rea{fu{w\x,t,e),x,e)} < -a 

uniformly for all (x, G E*, i e Z and small e > 0. 

H 3. In an internal layer, as ^ — > ±00 and e ^ 0, both and duf/d^ 
approach the corresponding values of w''^^ or w^^^ at common bound- 
aries. More precisely, if i = 2i, then for any > 0, there exist A^, eo > 
such that €q~^ > N, and for < e < eo, t > 0, 

\w\^,t,e)-w'-\f-\t),t,e)\<i^, for - e/^-i < ^ < -AT, 
\w'{t t, e) - w'+\y'{t),t, e) I < /X, for e^^-^ >^>N. 

Here the function w'' is expressed in the stretched variable C — ~ 

Let C = C(0 be a function of ^ such that 

'I, for < ef^-\ 

fore<-e^-S 
e^-^ for ^ > e^-^ 
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For each i > 0, i = 2£, consider the operator A'{t) : L'^{R) ^^(M), 

A\t)u = u^^ + Dtr]'{t, e)u^ + fu{w\^, t, e),7]'{t, e) + e|, e)^. 

H 4. A''{t), i = 2i, t > 0, has a simple eigenvalue A*(t,e) = eAg(t) + 
O(e^). The rest of the spectrum is contained in {ReA < — cx}, a as in 
H2. Moreover, for the limiting operator ^'(oo), we have, 

Ao(oo) < Aq < 0, uniformly for all i = 21. 

We now introduce a coordinate system 9^*(2) for the region E\ 

91(0)^91^(2) : E^^(x,^)^(e,T)eE^ 

Here the image of is denoted by S*. 

In the new variables the residual and jump errors are 

-g^ = w\-w\^-f\ (2.1) 

-nr) =(f.\ {w\[y' - y^-^]/2e, r, e) - w^+\y - y^+'\/2e, r, e)). 

(2.2) 

Here / = f{w\ \{y'~^+y')+^i^ e) with w' = w\i, r, e), and y^ = y^(er). 
The main result of this paper is the following 

Spatial Shadowing Lemma . Let {w*} he a formal approximation 
of solutions for (l.l). Assume that the Hypotheses H1-H4 are satis- 
fied. Then there exist positive constants Jq, J2, Sq and a negative con- 
stant 7 < — Cse satisfying the following properties. Assume that {w'} 
is a formal approximation near {w^}, with 

\w' - < Cie^\ i e Z, (2.3) 

for some Ci > and ji > 1. Assume that for the approximation {w*}, 
we have 

\9'\x(t,\-y) + l^1xO-75x0.25(R+^^) < C2e^\ j2 > jo- (2.4) 

ThenforO < e < eg, to any locally functionuQ with 1%— w*(0)|^i(2in{T 
€2^-^'^ , there exists a unique exact solution to (1.1) that satisfies Uexacti^, 0) 
Uq{x), and 

\Uexact ~ '*^*lx2.i(S%7)) — 0{e-'-^), i G Z. 

Here J3 = min{ji, j2 — J2} , J2 > is a constant that does not depends 
on e. All the norms in this lemma are expressed by the coordinate 
system 91^(2). 
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Remark We explain why the second approximation {w*} is needed 
in the lemma. Suppose that {w^} comes from a matched expansion 
of order e™. The constant jo is determined by Jo is stable with 

respect to small perturbations. Unfortunately it is not guaranteed that 
the residual and jump errors for {w*} are bounded by C2e^°. The 
formal approximation {w^}, on the other hand, can be obtained by 
adding higher order expansions to {w'^}- It will satisfy (2.3), therefore, 
share the same jo with {?!/}. Also its residual and jump errors, will be 
bounded by C2e^'^ for any desired j2 > if sufficiently many terms are 
added. 

Since we know how to compute formal approximations to any desired 
accuracy [25], the specific values of jo, J2 are not crucial. This was first 
observed in [4] . The constant ji is important since ja = ji if j2 is large. 
To compute ji, we can take {w^} as an expansion of order e"^+^, for 
higher order {w^} yields the same ji. 

Proof of the Spatial Shadowing Lemma. The proof is based on Theo- 
rems 5.1, 6.1 and 6.2. 

Let At, At, tf, r > he the constants as in Theorem 6.2. In the 
coordinates 9^(0), let Ij = [jAt, (j + l)At], < j < r - 1 and /oo = 
[tf, 00). A partition of the time axis 

j=0 

has been constructed in §6. Let 7 < be as in Theorem 5.1. 

At each t = jAt, < j < r — 1, or t = 00, with the X^'^ norm 
evaluated by the coordinates 9^*(2), 

\w'{i) - w\i)\^ < C\w' - w'\x2,i^j:\'y) < < Ce. 

Recall that the coefficients of the linear equations in §5 and §6 are 
obtained by linearizing at t. Thus the linearization around {w'^it)} 
and {w^it)} are e-close. The constant Ji in Theorem 6.2 depends on 
the coefficients of the linearization in a complicated way, but is stable 
with respect to small perturbations of the coefficients. Therefore, there 
exists Co > such that the same Ji is shared by all the {w*} satisfying 
(2.3). Let now jo = ^"1 + 71 + 1.5, where ri > 1.5. 

Since j2 > jo = ri + Ji + 1.5, all the terms in (2.4) are bounded 
by o(e^^"'"^'^+'^^). We now use the coordinates 9^*(1), as in §4, in each 
region E* H /j, < j < r — 1 or j = cxo. Recall that in the coordinates 
1>1*(1), n Ij maps to ^ x 1^ where = [0, At], and E* fl loo maps to 
Jl^ X R+. One can verify that the residual and jump errors are still of 
0(6^2) = o(e^i+"^i+^-^) in the coordinates ^^(1). After all, in the region 
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n Ij, the two coordinate changes 

9^(0) 

m{o) 9i'(2). 

only differ by a near identity mapping y = x — t,e). Therefore, 
in the finite time intervals Ij, Theorems 6.1 and 6.2 can be applied 
to obtain a unique exact solution, for each initial data uq satisfying 
\uo—w^{0)\jji(^Qi^^ < C2e^^. According to Theorem 6.2, the accumulation 
error on all these intervals are small so that the solution in the final 
interval I^d is guaranteed by Theorem 5.1. Denote the solutions by 
w^ + u^j, < j < r — 1, in n /j and + in E' n /oo. The solutions 
satisfy 

mm^i^.i.) <Cei-~'^-'-\ inQ} x7„ 

It is easy to verify that, with the left side evaluated by the coordi- 
nates ^\2) and the right side by IH^(l), 



r-l 



\u- 



X2,1(S%7) ^ Cel^l''^'^[^ |'u}|^2.i(Qixj^) + |M*oo|x2.i(Qi^xR+,7)]- 




Since r ^ log(l/(Coe^))/(e7Ar) and I7I < C^e, cf. §6, we have 

The number of the terms in the summation is r = o(e~^'^) if e is small. 
Thus, 

Let now J2 = Ji-|- Js+S. The proof of the lemma has been completed. 

□ 



3. BASIC Lemmas 

A function /(s) is in the Hardy- Lebesgue class ^^(7), 7 G R, if 

(i) /(s) is analytic in Re(s) > 7; 

(ii) {sup,>,(r^ \f{a + ^c.)prfa;)}V2 < 00. 

7^(7) is a Banach space with the norm defined by the left side of (ii). 
Based on the Paley- Wiener Theorem, [32], if e-'^f{t) e L'^{R+), then 
f{s) e 7^(7), vice versa. 

For k >0 and 7 e M, define a Banach space 

n'^ij) = {u{s) I u{s) and (s - j)''u{s) e 7^(7)}, 



\u 
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For any 7 e R, /c > 0, there exists C — C{j, k) such that 
C-\l + \s\^) < 1 + |s - 71'= < C(l + \s\''). 
Therefore an equivalent norm for 7i^(7) is 

/oo 
\u{a + + |(7 + iijj\^^)duj. 

_ . , -00 

It can be shown that if e-^'*/(t) G i/o^(M+), then /(s) e 7i^(7). 

Let n^^''^^{^) = :^'=i(7) X Ti'^^^^), Let i/J^"(R+) = //J(R+) x 
H'q{B+). Clearly, (e"^*?/, e"^*!-) G H^'^^W^) if and only if G 
'j-(_kixk2(^ryy Define a s-dependent norm in M^**. For {u,vy e M^**, 



= {l + \s\''')\u\ + {l + \s\''')\v\. 



An equivalent norm for {u,vy e 'H^'^^^'^i^) is 

/oo 
IKt')1|-x.,(,+,,)du;]^/'. (3.1) 
, -00 

Consider a linear equation 

ut = u^^ + Vi^u^ + ^(Oxi + /(t). (3.2) 

Here A{^) is a n x n matrix, continuous in ^, and T^(^) is a continuous 
scalar function in ^. Applying the Laplace transform to (3.2), we have 

su = % + V{Ou^ + A{^)u + />). (3.3) 
Convert it into a first order system 

Let T(^, C; s) be the principal matrix solution for (3.4) with s as a 
parameter. Let 5 be a subset of the complex plane and 11 C M be an 
interval. We say (3.4) has an exponential dichotomy in R^" for s G 5 
and ^ e f2, if for each s e 5, (3.4) has an exponential dichotomy in R^" 
for ^ e R. The projections PsH, s) +P„(,^, s) = / in R^" are analytic in 
s and continuous in ^. With the s-dependent constants K{s), a{s) > 0, 
we have 

C; s)Ps{C: s) = Psit s)r(e, C; «), C > C, 
m, C; s)Ps{C, s)|m- < i^(s)e-'^(^)i«-^i, e > c, 
m, C; s)Pu{C, s)|m- < K(s)e-'^(^)i«-^i, e < C- 

We say (3.4) has an exponential dichotomy in E^^^^'^{s) for s G 5 
and ^ e fi, if there exist projections Ps{i, s) + P„(^, s) = 7 in E^^^^'^{s), 
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analytic in s and continuous in ^, and constants K, a > 0, independent 
of s e 5, such that 

r(e, C; s)Ps{C, s) = Psit s)T{^, C;s), e > c, 

\t{^X;s)Ps{C,s)\e.,..,^,) < ;^e-"(v^+i)i«-^i, e > C, 
m,C;s)Pu{C,s)\EHx>^,^,^ < i^e-"(\^+i)i«-^i e < C- 

We say (3.4) has an exponential dichotomy in Ti^'^^^'^ {^) for ^ E Vt, 
if there exist projections Ps{C) + Pu{0 = in 7Y^^^^^(7), continuous in 
^, and constants K, a > 0, independent of s e 5. And if by specifying 
{u, v) ai ( E fl, solving (3.4) and denoting the solution map by T(^, (), 
we have 

(a) T(^, () : 7lPs{C) — ^ TZPg^C) is defined and continuous for ^ > C; 

(b) T(^, : 7lPu{C) ~^ T^PuiO is defined and continuous for ^ < C; 
|T(e,C)Ps(C)lH^ix.,(^)<Xe-"l«-«l e>C, 

|r(e,C)P«(C)k^ix^2w < xe-°i«-^i, e < C-' 

Lemma 3.1. Assume that (3.4) has an exponential dichotomy in £'0-75xo.25^g'^ 
for Re{s) > 7 and ^ e [a, b] . Then 

(1) (3.4) has an exponential dichotomy m 7i°''^^^°'^^(7) with the same 
projections derived from those in _£;0-75xo.25^_gj ^j^^ same constants 
K, a. 

(2) Assume that supJ|A(0| + \V{C)\] < M. Let (j) e 7i°-^5^°-25(7), 
{u,vy{i,s) = C-^{T{^,a;s)Ps{a,s)(t){s)). ThenuE H^'\[a,b]xR+,-f) 
and is a solution to (3.2) with / = 0. ^4/50 

Similar result also holds for {u, vY — JC~^{T{^, b; s)Pu{b, s)(f){s)). 

Proof Let e 7^0-^5x0.25 (^)_ Let w{^,s) = T{^X; s)Ps{C s^s), ^ > 
C > a with Ps{C, s) being the projection associate to the exponential 
dichotomy in the space ^0-75x0.25^^^^ r^Yie function w is clearly analytic 
for Re(s) > 7 and satisfies (3.4), with / = 0. 

Let s = a + iuj. To show (1), for each ex > 7, ^ > C, using (3.1), 

/oo p 00 

k(e,s)||o.r5xo.25(,)da; < / X2g-2a(l+|sr-5)(€-C)|^||^^^^^^^^^ 
-00 J —00 

<XV"(^-^)|0|^o.r5xo.25(^). 

This proves that w e 7y°-^^^°-^^(7) with the desired decay estimate. 
It is also clear that w e TIP si since the projection is derived from 
that in £!0-75xo.25^g^^ therefore, is commutative with the solution map. 
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This proves half of (1). The other half can be proved by considering 

w{t s) = T(e, C; s)Pu{C, s^is), ^<C<b. 

To show (2), let {u, vY = w{^, s) = T(^, a; s)Ps{a, s)0(s), ^ > a, and 
s = a + iu. For a fixed cr > 7, we want to show |'u(^, s)\ and \su{^, s)\ 
are in I/^([a, 6] x M), as functions of (^,0;). 

/oo rh roo rb 

/ {l + \s\^)\u\^dC(L;<C / / 
-00 Ja J —00 Ja 



< CK^ 
< 



/OO n 
-00 J a 



1 + |s 

b 



0.5- 



W I £;0 . 75 X . 26 (s) (i(x; 

0.5^^-2a(l+|s|0■B)(^-a) 



£0.75x0.25^5) 



d^duj 



2a 1^ lw°-^'5^"-25(7)- 
By the inverse Fourier-Laplace transform, we have 

|e~'^*1i|L2([a,6]xM+) + |e~'^*1it|L2([a,6]xIR+) < C|0|-ho.76xo.25(^) . 

Letting a — > 7, we found that u, Uf G L'^{[a, b] x W^, 7). Similarly we 
can show that |f |i2(^) = \u^\l2(^^^ < C|0|>^o.75xo.25(^). From the equation 
(3.2) itself and sup^{|A| + \V\} < M, \u^^\l2(^^^ < C|0|>^o.75xo.25(...),). This 
completes the proof of (2). □ 

A consequence of Lemma 3.1, is that if (3.4) has an exponential 
dichotomy in £'0-75xo.25(^_g-) f^j, Rc{s) > 7 and ^ G [a, 6], then (3.2) has 
an exponential dichotomy in i/Q^'^^"'^^(7) for ^ e [a,b]- That is, for 



any ^0 ^ [a, b] and (m°, v^) G Hi 



0.75x0.25 




(7) , there exists a decomposition 



(M°,t'°) = (u^.v^) + (M^,f°) such that (3.2), with / = 0, can be solved 
in [^0,^] with {u{^o) , u^{^o)) = It can also be solved in [a,^o] 

with {u{^o) , u^{^o)) = iuu,v^). Moreover, for some K,a> 0, 



MO. 



u 



In the next lemma, wc study Wr = u^^, u{^, 0) = and its Laplace 
transform in some detail. Consider 



-su = 0, se 6(M) = {| arg(s)| < 27r/3, \s\ > M}, 



(3.5) 



where M > 1 is a constant. Converting into a first order system, we 
have 



— su. 



(3.6) 



The eigenvalues are A = ±^/s and the eigenspace for each A is n- 
dimensional. We choose the major branch of ^/s. Thus | a,Tg{^/s) \ < 
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7r/3, Re^i > ^J\s\/2. Eqn (3.6) has an exponential dichotomy in ^ G - 
Moreover, the decay rates are e~V^'^'''^. One can easily verify that 

^vj~2\ y/su + v J ' 

^■u\ _ ^ ( U — v/^/s\ 
yv) ~ 2 \-^/sU + v) ' 



Puis) 



where and are spectral projections associated to the eigenvalues 
—^/s and ^/s respectively. 

For any A; > 0, using the s-dependent norm in R^", we have 

Ps{s) (l) = \ + + \V~SU-rv\{l+\s\')) 

\ J B(fc+0.5)xfe(5) ^ 

<^(|«|(l+|s|^+°-5) + |i)|(l + |s|'=)). 

Here we have used the fact |s| > 1. Thus |Ps(s)|^(fc+o.6)xfc(s) < 1.5. 
Similarly we can show that |Pu(s)|^(fc+o.5)xfc(5) < 1.5. Observe that for 
|s| > 1, ^ > 0, we have e'^^/^ < g-(i+^)«/4_ ^Ve have thus proved 
the following 

Lemma 3.2. Eqn (3.6) has an exponential dichotomy in k > 

0, for s e 6(1) and C&R, with K = 1.5 and a = 0.25. 

2ti 

Remark The angle — in the definition of the sector can be replaced 

TT 

by any — < 6 < tt. However, the constants K and a will be different. 

We shall often use the well-known Roughness of Exponential Di- 
chotomy Theorem. See Coppel [6] for a proof for / = M-|-, M~ or M. 
The proof extends to the case / being a finite interval. Below, we 
generalize the lemma to the space £;('=+o-5)x*:(s)^ A; > 0. 

Lemma 3.3. Assume that an ODE in R^'^, 

D^Y = A{C,s)Y 

has an exponential dichotomy in E'^^~^'^'^^^''(s), k > for ^ E I , and 
s E S, where S G C, and I can be M, M"'", M.~ , or a finite interval 
[a, b]. Let the projections be Ps{C, s) + Pu{C, s) — I, and the constants 

be K, a > 0. Then there exists cq > such that if S ^ sup^^i\B{^)\ < 
Coa(l + M°-5) for some M >0, then 

D^Y^[A{t)+B{t)]Y, 
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also has an exponential dichotomy in with ^ E I and 

s e <S n {|s| > M}. Let the projections for the new system be Ps{C, s) + 

PuiC^ s) = I and the constants K, a > 0. Then there exist Ci, C2, C3 > 
such that |P,(^,s)-P«(^,s)|£;(fc+o.B)xfc(^) < Ci5/(a(l + M^-^)), |Q;-a| < 
C25/{1 + M^-^), andk < C3. 

Proof. Following the idea of [6], or [26], Theorem 3.5, we first obtain 
projections to the stable and unstable spaces for the perturbed system 
in £;('=+o-5)x'=(s). Solutions on these spaces decay like e-("(i+l^l° ')-^2<5)|?-ci^ 
if the initial point is C- However, 

a{l + \sr) - c,5 > a{l + |.|°-^) - ^^^(1 + l^H- 

C2S 

This explains that the decay rate has the form {a — - — + 

\s\'-') ^ □ 

Lemma 3.4. Consider (3.4) with sup^{\A{()\ + \V{C,)\} < a < 00. Let 
k > be a constant. Then there exists M > 1, depending on a, such 
that (3.4) has an exponential dichotomy in £;(*^+o-5)x*;(s) for s e &{M) 
and ^ e M. 

Proof. From Lemma 3.2, system (3.6) has an exponential dichotomy 
with K = 1.5 and a = 0.25 in the region s G 6(M), M > 1. We then 
use Lemma 3.3. We can choose M larger so that S — a < 0.25co(l + 
M°-^). The result follows from Lemma 3.3. □ 

Corollary 3.5. (1) Consider (3.4) with the same hypotheses as in 
Lemma 3.4- Let S G C be a closed subset. Assume that iS fl {|s| > 
M} C (3(M) for some M > 0. Assume that (3.4) has an exponential 
dichotomy in M^" for s E S and ^ eM.. Then (3.4) has an exponential 
dichotomy in E^''+'^-^^'"'{s), k > 0, for s e S and ^ & 

(2) Let (Jq G M. Assume that (3.4) has an exponential dichotomy in 
M^" for Re{s) > ctq and ^ G M. Then it has an exponential dichotomy 
in ^(*^+o-5)xfc(s), A; > 0, for Re{s) > ctq and ^ G M. 

Proof. Let M > 1 be the constant as in Lemma 3.4. We can choose M 
larger so that «S n > M} C 6(M). By Lemma 3.4, (3.4) has an 
exponential dichotomy in £!('=+o-5)xfe(s) for s G Sr)&{M) with constants 

Ki, ai. 

In the compact subset 5 fl {|s| < M}, (3.4) has an exponential 
dichotomy in M.^"" with constants K2, 0.2- The norm of E^^^'^'^^^^i^s) is 
equivalent to the norm of R^" uniformly with respect to s, \s\ < M. 
Besides, 

X^e-^^l^l < X2e-"^(^°"+^)-'(l^l°"+^)l^l, if \s\ < M. 
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Therefore, we find that (3.4) has an exponential dichotomy in £'('=+o-5)x'=(s) 
for5n{|s| <M} and^eM. 

By choosing larger K and smaller a, the result of (1) follows. 

Observe that for any fixed (Tq, if M is sufficiently large, we have 
{Re(s) > (Jo} C 6(M). The result of (2) then follows from that of 
(1). □ 

Remcirk Corollary 3.5 is valid only the exponential dichotomies are 
considered for ^ e R. Since it uses the "uniqueness of exponential 
dichotomy for ^ G M". Caution must be exercised when using it to 
prove, for example. Lemma 4.1 later. 

Lemma 3.6. Assume that A is a constant matrix and V is a constant 
scalar. Suppose that there exists (Tq > such that Rea{A) < — (Tq. Then 
for any < S < ao and k > 0, (3.4) has an exponential dichotomy in 
^(fc+o.5)xfc(g) /or Re{s) > -ao + S and ^ e R. 



Proof. Let A be an eigenvalue for J(s) = ^ yjj . Let x + iy 

be an eigenvalue for si - A. Then A = {-V ± ^/V'^~+~^x+~^^y)/2. 
Using the fact Re(s) > — ctq + 5, we find that counting the multiplicity, 
there are n eigenvalues with 

Re(A) < {-V - VV^ + A5)/2 < 0, 

and n eigenvalues with 

Re(A) > {-V + VV^ + A5)/2 > 0. 

Therefore, (3.4) has an exponential dichotomy in R^'* for Re(s) > —ao+ 
S and ^ G M. From Corollary 3.5, it has an exponential dichotomy in 
^{k+o.5)xk(^g^ for the same s and ^. □ 

Consider 

= [sI-A{eO -m)/) © + (/)• ^^-^^ 
When e > is small, (3.7) is a system with slowly varying coefficients. 

Lemma 3.7. (Existence of exponential dichotomy in regular layers) 
Assume that the matrix A{^) and the scalar V{^) are hounded func- 
tions for ^ e R. Assume that 

and there exist ctq > such that for all ^ G R; 

Re{a{A{i))} < -ao < 0. 
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Let < 5 < (Tq. Then there exists eo > such that for < e < Cq, 
there (3.7) has an exponential dichotomy in k > 0, for 

Re{s) > -(7o + S and^ eR. 

Proof. According to Corollary 3.5, (2), we only need to show that (3.7) 
has an exponential dichotomy in R^" for any s e {Re(s) > — ctq + S} 
and ^ e M. 

Fixed such (s, ^) as a parameter, from Lemma 3.6, system (3.7) has 
an exponential dichotomy in R^'* as a system with constant coefficients. 
Because of the bound on \A\ + \ V\, the constants K, a can be chosen 
independent of ^, see [6]. We then use Proposition 1 from [6], page 50. 
Since 

\D^A{eO\ + \D^V{eO\<Ce, 

if eo > is sufficiently small and < e < eo, all the conditions in that 
proposition are satisfied. Therefore, (3.7) has an exponential dichotomy 
in for s e {Re(s) > -ao + 5} and ^eR. □ 

Lemma 3.8. Consider a closed liner operator A : L^(]R) — > L^(R) 
defined by 

Au = u^^ + V{i)u^ + A{i)u. 

Assume that (3.4 ) has an exponential dichotomy in £''^'^^°(s) for Re{s) > 
—(To and ^ e R; then A is a sectorial operator with 

l(A-^)-^|< ^ 



l + IAI' 



for all ReX > — cq. 

Proof. For g e I/^(R), Re A > — Uo, we want to solve 

Au — \u — g. 
Convert it into a first order system, 

0^ = (a/ -A{i) -v\i)l) + • 

From the existence of exponential dichotomy, the only solution can 
be written as 
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Since e ^0-5x0 (A), we have 

<i^e-(|Ar+^)l«-CI|^(^)|^„, ^>^. 

Similarly, 



We now have 



U \ , , I X i0.5\-l 



|m(OIiR" < I y^j |£;o.5xo(A)(l + |A 

< (1 + |A|°-^)-^( + r);re-(^-^i^i°-^)i«^^i|^(c)|dC 

J — oo J 5 

< X(l + |A|0-S)-i[e-"(i+W°")l«l * |5(0|]. 

Here "*" denotes the convolution in ^. Prom a standard estimate on 
convolution, 

Ix^U. < X(l + |A|°-^)-^|e-(^+W°-^)l«IUi|^(0|L^ 
< 2X(l + |A|°-^)-2o;-^|^|i2 
<2X(1 + |A|)-V-^|^|l2. 
The last estimate indicates that A is sectorial in (R) . □ 



Lemma 3.9. LetS dChea closed subset andSn{\s\ > M} C 6(M) 
/or some M > 0. Assume that (3.4) has an exponential dichotomy in 

^2n s E S and ^ G M~ or M"*" respectively. Let an operator A be 
defined as in Lemma 3.8. Assume that S C p{A). Then (3.4) has an 
exponential dichotomy in ij('^+0-5)x'=(s) for s E S and G M. 

Proof. From Corollary 3.5, one only has to prove that (3.4) has an ex- 
ponential dichotomy in M^"^. Since the exponential dichotomies exist in 
and respectively, one only needs to show the transverse intersec- 
tion of TZPs{0~^, s) and TZPu{0~, s). If for some sq £ ^, this is not true. 
Then there exists 7^ such that 4> G nPs{0+, sq) n nPu{0~, sq). Let 

= r(^,O;so)0- Then u = C'^u G L'^{R) and satisfies SqU = Au. 

This is a contradiction to Sq G p{A) . □ 

Lemma 3.10. Assume that 7 G M, and (3.4) has an exponential di- 
chotomy in R^" for Re{s) > 7 and ^ G R" or R"*" respectively. Assume 
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that A (see Lemma 3.8) has a simple eigenvalue Xq, Re{\o) > 7, and 
no other eigenvalue in Re{s) > 7. Then (3.4) has an exponential di- 
chotomy in ^(fc+0-5)xfe(^^5) {Re(s) > 7} n {\s - Ao| > 5} and ^ G M. 
Here 6 is any positive constant. Let the constants of the exponential 
dichotomy be K^, as. Then as can be independent of 6, but Ks ^ 00 
as 5 — > 0. More precisely, 

7eP,(0+, s) © 7^P„(0-, s) = M^", (3.8) 

for all {Re{s) > 7 and s ^ Xq}. The projections Pi + P2 = I, defined 
by the above splitting have a pole of order one at s = Xq, i.e., 

\Pj\ < C(l + \s - Aol-^), J = 1, 2 for Re(s) > 7. 

Proof. We first show that {Re(s) > 7 and s 7^ Aq} C p{A). Let s 
be a point from the left hand side. Then (3.8) is vahd, or s e cr{A), 
which is not true since s 7^ Aq. Prom (3.8), (3.4) has an exponential 
dichotomy for ^ G R. Using the same integral formula as in Lemma 3.7, 
we find that for each g E L^(M), there exists a, u E i7^(]R) such that 
Au — su = g. Thus, s G p{A). 

The closed subset {Re(s) >7n|s — Ao|>5}c p{A). Thus from 
Lemma 3.9, (3.4) has an exponential dichotomy in E^^^^-^'>^^ [s) for 
s G {Re(s) > 7 n |s - Aol > 5} and ^ G M. 

Consider |s — Ao| < 5, 5 small. Based on a lemma due to Gardner & 
Jones [19], There exist n independent vectors {w^{s)}'^^^ that form a 
basis for TZPu{Q~,s) and n independent vectors {wf{s)}'^^^ that form 
a basis for TZPs{Q'^, s). The vectors are analytic in s, so is 

D{s) = dct{w-{s) . ..w-{s) w+(s) . . .w+(s)). 

The results from [19] also assert that 

(1) D{s) = in |s — Aol < 5 if and only if s G crp(^); 

(2) The order of the roots ol D X — Xq equals the algebraic 
multiplicity of Aq as an eigenvalue of A. 

We infer that s = Aq is a simple root of D since Aq is a simple 
eigenvalue. 

For w G R^", |w| = 1, consider 



w = {wi (s) . . . {s) wi{s) . . . w+(s)) 



\C2n{s)) 



From Cramer's rule, we easily find that Cj{s) = 0{\D{s)\ ^) = 0{\s — 
Ao|~^). Therefore, PiW = J2i=iW^{s)ci{s) = 0(|s — Ao|~^). Similarly 
P2W^O{\s-Xo\-^). 
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The desired estimates on Pi and P2 follow by combine estimates in 
{Re(s) > 7 n |s - Ao| > 5} and {|s - Ao| < 5}. □ 

Lemma 3.11. Assume that A is the operator defined in Lemma 3.8 
and is sectorial in L^(]R) with constants ckq > 0, C > such that 

S = {| arg(A + ao)| < 0, n/2<e< 7r}\{-Q;o} C p{A), 
\{X-A)-%2<Co\\ + ao\-\ forXei:. 

Let 7 < with < ao/4. Let g e L'^{R x R+,7) and Uq G H'^{R). 
Consider the initial value problem, 

Ur = Au + 

«(e,0) = Mo(0, eeM, r>0. 

Then the solution u G H'^'^iR x IR+, 7), and 

< C(|iio|i/i + |£/|l2(^)), 

where C = a/2Co(1 + ^- ). 

cto + 37 

Proof. B — A — ^ is a, sectorial operator with 
\{X-B)-%2 < 



Co 



|A + ao + 7I 

for A e (I arg(A + ao + 7)| < 6*, 7r/2 < 9 < 7r}\{-Q;o - 7}- Suppose 
now Re(A) > 7, 

|A + ao + < [(A - if + K + 27)T'/' 

<^[|A-7| + («o + 27)]-' 

< V2[|A| + ao + 37]"' 

<(l + (ao + 37)-') ^ 



Al + l 



Therefore 



|(A-S)-^U2<C^-^, (3.9) 

where C is the constant defined in the lemma. 
The solution of the initial value problem satisfies 

e-''*u^e^*UQ + e^* *G{t), 

where "*" denotes the convolution in t and G{t) — e~'^g{t), with 
1^11,2 = |^|l2(^). Prom (3.9), and [28] 

|e^**G(t)|if2,i < C\G\l2, 

|e^*^io|i^2.l < C\uo\m- 
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The desired result then follows easily. □ 

4. Coordinates Change and Linearization 
We first construct new coordinates in the region 

n Ai], i > 0, or E^n[i/,oo), 

to straighten the boundaries F^-^&P. Here At — eAr, At > and 

tf > 0. Let < 9{x) < 1 be a function, ^(0) = 1/2, e{x) = 1 for 
x > 1 and e{x) = for x < -1. Let m\t) = {y'-\t) + y\t))/2 and 
Ay\t) = y\t)-y'-\t). For each fixed t > or t = oo, let / = [t, t+At] 
or I = [tf, oo) respectively. Define a change of coordinates in fl /, 

y = y*{x,t,t,i) = X -E{x,t,i,i), 

Note that y^{t) really depends on e, so is y* and S. For typing conve- 
nience, we will also drop t.i iny*, S if no confusion should arise. 

From its definition, it is clear that S = for t = t. 

Let C > be a given constant. For the case i = oo, from §2, HI, 
there exists tf sufficiently large such that for / = [t/, oo), 

\y\t)-y\t}\ + \y'-\t)-y'-\t}\<C, tel, (4.1) 

For the case < t < oo, if e > is sufficiently small, then for / = [t,t + 
At], we also have (4.1). If C is sufficiently small, with the corresponding 
tf or At, we have 

d^E{x,t,t,i) <^, {x,t)eE'r\I. 
Thus, y — y*{x,t,t,i) has a smooth inverse x — x*{y,t,t,i), y^~^{i) < 

y < fit), t e I. 

The change of coordinates has the following properties: 

{1) y = X a.t t = i; 

(2) A neighborhood of F* in E' fl / is shifted by the amount y^{t) — 
y^(t)- Thus in the new coordinates, F* = {y = y^{t)}. For the same 
reason, F^^^ = {y = y''\i)}, and S^n / = {t/'^^(f) <y< y\i)}. Also, 
if \x - y\t)\ or \x - y''\t)\ < Ay\t)/A, then 

a^S(a;,t,t,i) = 0, k>l. 

(3) The curve that bisects the region is straightened, y*{m^{t),t) — 
rrf{t). 

(4) (a) In a finite interval / = [t,t+At], \d^d^x*{y,t)\ + \d^d^E{x,t)\ < 
Ckh uniformly with respect to t > 0. 
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(b) In the interval /oo = [tf,'^), \9ydtX*{y,t)\ + \d^dtE{x,t)\ < 

The proof of (4) (a) is based on the continuity of D^y'^{t). The proof 
of (4) (b) is based on \Dy''(t)\ < Ce~^*, which in turn is based on HI. 

We now choose tf larger, if necessary, so that there exists an integer 
r > such that tf = rAt. Let Ij = [jAt, {j + I) At], < j < r - 1 and 
-^oo = [tf, oo). A partition of the time axis is made by /j. < j < r — 1 
and /oo- Define a change of coordinates hjy = y*{x, t, t, i) with i = jAt 
in n Ij and t = oo in E* fl /oo- Using the new variable, a solution of 
(1.1) in regular or singular regions satisfies 

u{x,t) = U{y,t) = U{x -E{x,t),t), 

e(Ut - UyEt) = e2((l - E^fUyy - E,,Uy) + f(U, X*, e). 

We now construct a new coordinate system in the region fl/j, < 
J < r — 1 or J = oo. Let t — jAt in E* fl Ij, including j — oo. Let 
Lj(e) = Afit}/2e and = (-L}(e), L^(e)). 

^(0) ^ : E^ n 3 {x,t) ^ (e,r) G Q} x 7^, 

y^x- E{x,t,t,i), 
e = b - m'it)]/e, 

f {t-jAt)/e, jy^oo, 

, [0,At], inE^n7,-,0<i<r-l, 
\ [0,oo) inE^n/oo- 

Using the coordinates we denote the solution in fi^ x < 

j < r — 1, or oo by r). The equation for W is 

[/; = + + ee,t),t)f/i + /([/\a;*K(t) + ee,t),e) 

(4.2) 

+A/?(c/\e,t). 

Here A/^ = if i = f. A/^i* is hnear in U': 

Ml{U\ i, t) = [-2E,{x*, t) + El{x\t)]Ul^ - eEUx\t)Ul 

+ [Et{x\m\i) + ti,t),t) -Et{x\m\t) + ei,t),i)]Ul. 

Here x* denotes x* (m* (f) + e^, t) and ^ = f + er if f = j At or t = ^/ + er 
if t = oo. 

Using (4), (a), it is not hard to verify that if t = jAt, < j < r — 1, 
then 
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Using (4), (b), ii t — oo, j — oo, then for any 7 < 0, 

Let t/* = + M*. Linearizing (4.2) around w'^{m^{t) + e^,t, e), where 
is the formal approximation, we have 

(4.3) 

Here A\i) = fu{w\m\i) + e^, t, e), m^(t)_+ ee,e), ^^(0 = Si(m^(t) + 
e^,t). One can verify that V'^^{^) — r]f{t,e) is the speed of the wave 
front, and is independent of ^. l^^^+i depends slowly on ^. — if 

i = 00. Denote uf{t) = tv\x* {m'- {i) + e^, t),t, e), 

A/-^ = jV;' + + [/(^/;Xt) + ti\ x*(m^(i) + eC, i), e) - /(«;^(i), x*(m^(i) + e^, t), e) 
-fu{w\t),x*{m\t} + eU),ey] 

+ [Uw'it),x*im'{t) + eU),e) - fu{w\t),m\t} + eC,t},e)]u\ 
When < j < r — 1, we have, in ^ x 

\Af%2<\g%. + C{At\u%., + \u%,.,), 
\J\f'{u') < C{At+\u%2.i + |m'|h2,i)|m' -mV2,i.(4_4) 

When j — 00, we have, in x 7^, 7 < 0, 

|AA^(«') - AA^(«)U(7) < C{e-^'f + |Mix2.i(7) + I^U2.i(^))|m^ - ^ix^.^^) 

We are led to solve system (4.3) with t = j At, = [0, At] if < 
j < r — 1, or 7:^ = M"*" if j = 00. The nonlinear term A/"* satisfies (4.4) 
or (4.5). The solution satisfies jump conditions (1.8) in the form 

(ij [Ame),r,e)-u^+\-L^\e),r,e)] = S\r), (4.6) 
initial conditions 

«'(e,0) = «^(O, (4.7) 
and compatibility conditions 

7ri<5XO) = vUme), 0) - xi^+^(-Lf ^(e), 0). 

Here tti : R^" — M", (m, i')^ — > m, for u,v ^ M", is a projection. For 
5* = (^1,5^)^ e f/'0-75xo.25^jj^^ observe only the trace of d\ is defined 
at T = 0. Therefore, the compatibility condition is posed only on 
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Let US discuss properties of (4.3), in regular or internal layers, and 
in finite or infinite time intervals. Assume the the coefficients of (4.3) 
have been extended continuously from ^ e O*- to ^ e M by constants. 
This is equivalent to introducing ^ in internal layers, see §2. If, af- 
ter freezing the coefficients /1*(C) and V*(^) at ^ = (4-8') below 
has an exponential dichotomy for ^ G M, then the projections for the 
dichotomy are denoted by -P](^o) and P^{C^o). 

Let A^u = u^^ + V^{C,)u^ + A^[^)u. The dual system associate to 
Ur — A^'u is (4.8'). When i = 2£ — 1, for each fixed A*({) satisfies 
H2, and \d^A\i)\ + \d^V\i) \ = 0(e). Prom Lemma 3.7, we have 

P 1. In the regular layers, (4.8') has an exponential dichotomy in 
£;(fe+o.5)xfc(^^^ for Re(s) > -a + cti, cti > 0, and C e M. Also the 
projections of the dichotomy at any £ 1^ are close to those obtained 
by freezing the coefficients, i.e., P;(^o, s) = P;(Co, s) + 0(e), P^(^o, s) = 
Pl(i,,s) + 0{e). 

Based on Lemma 3.8, we have 

P 2. In the regular layer, |(A - ^T^U^ < C/(l + |A|), for Re(A) > 
—a + (Ti. 

For i = 2i, if the coefficients are constants, A^"^^ {— L'^j^^ (e)) and 
y'+^(— L*"'"^(e)), system (4.8') has an exponential dichotomy for G M, 
due to H2 and Lemma 3.6. Similar assertions holds for (4.8') with 
coefficients A'-\Lf-\e)) and V'-\L^\e)). Using Lemma 3.3, and 
H3, we can show 

P 3. In the internal layers, there exists a large constant N > such 
that the dual system (4.8') has an exponential dichotomy in i^('^+0-5)x'^(s) 
for Re(s) > —a + ai, ai > 0, and ^ G [— oo, —N] or [N, oo) respectively. 
The projections are close to Pj~^(L*~^(e), s) and P^~^(L*~^(e), s) for 
^ < -N. They are close to P^+^-L^^e), s) and Pl+\- L'^^e) , s) 
for ^>N. 

Based on P3, (4.8') has exponential dichotomies in R^"^ for ^ G M~ 
and respectively. They are natural extensions form the dichotomies 
in (— oo, —N] and [N, oo). Prom H4 and Lemma 3.10, we have 

P 4. In the internal layers, (4.8') has an exponential dichotomy in 

^(fc+o.5)xA-^^.) ^ ^ >_a + 5n{\s-X'{e)\>S}} for all 5 > 

and ^ G M. The projections defined by the splitting 



7^P](o+, s) e 7^p^(o-, s) = r- 
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are bounded by C{1 + - ). Further more A*(e) = e\l,{i)+0{e^). 

In the hmiting case, we also have Aq(oo) < A < 0, i G Z. 

Remark Exponential dichotomies arc not unique in semi-infinite inter- 
vals [A^, oo) and (— oo,— A^]. In particular, the projections that define 
the exponential dichotomies in P3 and P4 are difi^erent. However, the 
unstable subspace at —N, 71P^{—N, s), and the stable subspace at N, 
TZPl{N, s), are unique. 

We now present an important lemma that is used in §5 and §6. 
Consider a sequence of equations, i £ Z, 

< = 4^ + V^'(04 + ^'(eK, ee^.re/^ (4.8) 

1.^(^,0) = 0, (4.9) 



(4.10) 

The compatibility 7ri5*(0) = is clearly satisfied. 
The dual systems for (4.8)-(4.10) are 



(4.10') 



where 5' e 7^0.75x0.25 Let 



B^sup{\A\0\ + \V\0\}- 

We first derive a Gronwall type inequality that governs the growth 
of solutions of (4.8'). A solution (j){^,s) of (4.8') satisfies the integral 
equation, 

0(e, .) = e'^^HiO, s) + e^(^)(^-f) (^_^°^^^ -v\c)l) 
Consider s G 6(M), M = 1 first. See (3.5) for the definition of 6(M). 
It is known that J{s) — has eigenvalues A = ±-\/i, each has an 

n-dimensional eigenspace. The spectral projections to the cigenspaces 
are bounded uniformly with respect to s in the norm of E^''^^^^''^^[s). 
Thus for ^ > 0, we have, 

m,s)\<Ce^\^\<t>{^,s)\+C f e^(^-^)S|0(C,s)|dC, 

^0 
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where the norms are in E^-'^^^^-'^^{s). The constants C are related to 
the norms of the spectral projections. Prom Gronwall's inequality, 

m,s)\ < Cie(v^+^^)«|0(O,s)|, e > 0, £>! = CB. 

(4.11) 

Now consider s e {Re(s) > — a-}\6(M), a compact set in C. By 
Gronwall again, we have, in the R^" norm, 

m,s)\<C2e''-^<t>(0,s)\, e>0. (4.12) 

The constants Ci, C2, -Di, -D2 are all positive. Results similar to (4.11) 
and (4.12) can also be derived for ^ < 0. By combining (4.11) and 
(4.12), we have, 

(4.13) 

for ^ G M and Re(s) > —a. Here the constants C and rj depend on a 
and B only. 

Lemma 4.1. Assume that (4-8) satisfies PI and P2 in regular layers 
(i = 21—1), and satisfies P3 and P4 in internal layers (i — 21). 
Assume either one of the following conditions is valid. 

(1) J = oo,Ii = M+, 7 = 7oe, 70 < anc? 7 - sup,{|A^(e)|} > de 
for some Oi > 0. 

(2) < j < r — 1, — [0, At]; 7 is a constant, independent of e, 
and —a < 7 < 0. 

Then, in any of the two cases, there exists eo > such that for < 
e < €0, (4-8)-(4-10) has a unique solution {u^}'^^ with e if^'^(f2*- x 
li, 7). Moreover, 

Proof. Consider case (1) first. Let — cto = —a + ai < inf{A*(e)}. Then 
—(To < 7. By PI, (4.8') has exponential dichotomies in regular layers 
for Re(s) > -ctq and ^ G M. 

By P3, in internal layers, there exists exponential dichotomies for 
Re(s) > —(To and ^ G (—00, — A^] or [N, 00). The range of the projec- 
tions at the boundary T\ nPs{-Lf^{e), s) and 7^P^(Lj(e), s), are close 
to the stable and unstable eigenspaces of a constant system. Thus, the 
projections defined by the following splitting 

nPl-^\-Lf\e),s)®nP:{L){e),s), 

are uniformly bounded with respect to e. For 5' G 7^0-75x0.25^^^^ 
there exist unique 0t(L^(e), s) G nPi{]Jj{e), s) and 0^+^(-L}+^(e), s) G 
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■)IhO-75x0.25(^) + |0*+^(-L}+^(e), •)|^0.75X0.26(^) < C | 5' | ^0.75x0.26 (^) . 

In a regular layer, let 

^ s)<f>i{-Liie), s)+T% me), s)0t(L}(e), s), \C\ < ^(e). 

In an internal layer, first consider 

(/. = r(e,-4(e),s)0:(-L}(e),s), ^e(-Li{e),0]. 

Since (4.8') has an exponential dichotomy for ^ G (—00, —N] or [iV, 00), 
therefore u — >C~^(7ri0), ^ e (— L*(e), — A?") is a solution of (4.8) with 
u e H'^'^{[—Lj{€), —N] X /^,7), according to Lemma 3.1. Moreover, 

|0(-iV,s)Uo.r5xo.25(,) <Xe-'^(^^+^)(^5W-^)|0;(-L}(e),s)Uo.75xo.25^ 

We want to show u e H'^'^{[—N,0] x /^,7) also. The exponential 
dichotomy used above has not been extended to [— A'^, 0]. However, 
using inequality (4.13), we can show 

|0(^, S)|£;().75X0.25(5) < C'e''^V^+^^^^+^)|0( — A^, S) lijO. 75x0. 25(5-) , —N < 1^ < 0. 

Let s = a + iuj with cr > 7. Let e be sufficiently small such that 
a{Lj{e) — N) > rjN. We have, as a function of u; and ^, 

/ {\s\^ + l)\u{ts)\''d^dw 

-00 ./-AT 

/oo ("0 
/ (|s|°-^ + l)|0(e,s)||o.75xO.25(,)det^C^ 

- II + «)||o.75xc.25(,)rfec^a; 

<C r e2''(\^+^)^X2e-2"(\^+^)(^^(^)-^)|0:(-L5(6),.)||o.75xo.25(,)rfu; 

/oo 
|0^(-L}(e),s)||o.75xo.25(,)O?a; 
-00 

•^0.72x0.25 

Here Co < e2''^-2"(^}(^)-^) ^ as L^e) ^ 00. 

Let[/He,5) =T^(e,-i^K^)'^)<^s(-^;(e)'«)' U\^,s) = T\^,me),s)(t>i{Li{e),s). 
We have shown that C'^niU^ e H^'\[-Li{e),0] x 1^,7). Similarly, 
C-^TTiW^ e i/2,i([o,L;.(e)] X 7:^,7). From P4, when Re(s) > 7 > 
supjA'(e), there exists an exponential dichotomy for ^ e M. This is 
not the same exponential dichotomy guaranteed by P3, see the remark 
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after P4. Denote the projections associated to that dichotomy by Pg 
and Pu- Let 

0:(O+, s) = P:(0, s)iU\0, s) - U'iO, s)), 
<P:{0-,s) = P^{0,s){U\0,s) - U\0,s)). 
We have, in i^o.75xo.25^_g^ norm, 

i0:(o+,5)i + ic(o-,^)i 

< C(l + L^)e^^-m(^)-^)\{S% 

\s - A'(e)| 

Consider = 0, s)0^(O+, s), ^ e [0,iv}(e)]. Using the same 

exponential dichotomy in [0, L^(e)], guaranteed by P4, 



m,s)\ < ;re-"(VH+i)«|0;(O+,s)|, Re{s) > 7. 
From Lemma 3.1, u = C^^7ii[Pi(p] is a solution in H'^'^{'j). 



\u 



/fM(^)<ci0:(o+,s) 

< (7(1 + L^)e^^-"(^5W-^)|{5^}|j 



-— ' I t« 1 1^^0.75X0.25 (^).^4_^g^ 

Similar result holds for = r^(^, 0, 8)0^(0", x), ^ e [-Lj(e),0]. 
For Res = (T > 7, let 



0^ 



r'(e,-L5(e),s)0:(-L;.(e),s)+T^(e,O,s)0UO-,5), ^ e [-L} (e) , 0] , 
0, s)0:(O+, s) + T^(e, L}(e), s), e e [0, L}(6)]. 

The function -u* = £^^7ri[Pi0*] is a solution for (4.8), but the jump 
condition (4.10) is not satisfied. In fact, (4.10') has an error 

|r (L}(6), 0, s)0:(O+, s) - T\-Lf\e), 0, s)(t^^\Q- , s)|^o.r5xo.25(,) 
< C(l + ;^-i^)e(''+'^)^-^"^5(^)|{5^}|^o..5xo..5(,). 

Due to the fact 7 — |A*(e)| > Oie and i^j(e) = e^~^ in internal layers, 
the error is much smaller than {5*} if e is small. It is clear that the 
desired solution that satisfies both (4.8) and (4.10) can be obtained 
by iterations. The initial condition (4.9) can be verified by the Paley- 
Wiener Theorem. The estimates on |-u*|j:/2,i('y) also follow easily. 

We now consider case (2), < j < r — 1. The beginning part of 
the proof is identical to that of case (1). since it is possible that 7 < 
supi(A*(e)), We do not have (4.15) for Re(s) > 7. Since A*(e) = 0(e), 
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if e is small, (4.15) is valid for Re(s) > 1. It turns out that (4.16) can 
be replaced by 



\u 



< C(l+ ^^^^)e''^-"(^M^)-^)|{57|^o.r5xo..5(,). 



Define 0' as in case (1). We need a more precise estimate for u\ In an 
internal layer, when ^ e [— L^(e), 0], we can write it' = u\ + u\. Here 

< = £-^[7ri(r (e, -L;.(e), s)0:(-L}(6), s))l 
4-£-i[7ri(r(e,0,s)C(0-,s))]- 

Similar to case (1), |it^|^2,i(-y) < C|{(5*}|ifo.76xo.25(^). If e is small such 
that 1 - |A^(e)| > 1/2, we have 

< Ce^^-"(^^(^)-^)|{5'^}|^».7sxo..5(,). 

Using the fact 7:^ = [0, Ar] is a finite interval of length At, we have 

Ul* I . . < p(l-7)Ar I i I 

l"'2li?2,l(Q»x7^,7) — I"'2lif2,l(n»x4,l)- 

We have shown 

By the Trace Theorem, the boundary value of U2 ~-^*(^) satisfies 

|4(-L}(e), •) 1^,0.76X0.25(4,^) < Cie(^-^)^^+(''+")^-'^^5(^)|{(5^}|^0.76x0.25(^), 

Let < e < eo be small such that 

The smallness of eo depends on At and 7. We then have 

|«2|h2.1(7) - C'|{'^*}|//0-75x0.25(^). 

l'"2(~-^*(^)) ■)|j/0-'^5x0.26(^) < -|{5'}|if0.75x0.25(^). 

Similar result can also be derived for = u\ + in [0,L*(e)]. Of 
course, all the above are trivially valid in regular layers. 

The rest of the proof is just like that of case (1). -u* = £^^[7ri0*] is 
a solution to (4.8). The jump condition is not satisfied, but the error 
is small. The precise solution {u^^°^^ can be obtained by successive 
iterations. □ 
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5. Corrections in the Final Interval [t/,oo) 

As shown in §4, the original problem is equivalent to the nonlinear 
system (4.3), (4.6) and (4.7). In this section we study this system in 
the infinite interval I^o — [^/^oo)- Wc have shown that (4.3) satisfies 
PI and P2 in regular layers, and satisfies P3 and P4 in internal layers. 
The nonlinear term A/"* satisfies (4.5). Our main result is the following 
existence theorem which also contains an estimate on the solution in 
terms of the initial, residual and jump errors. 

Let 7 = sup{— €7, — for all i e Z} where 7 is the constant in 

(4.5). Let e X{Q^^ x M+,7), u], G H\D^^), r = {t - tf)/e and 
e xo-75xo-25(7). 

Theorem 5.1. Letri > 1.5, \g'\x{-y) = o(e''i+^-^), |Molffi = o(e''i+°-^), |5*|jyo.75xo.25(^) = 
o(e''^). Let tf be e dependent such that e~'^^f < CqC^ , where 7 is the 
constant in (4-5). Then there exists Cq > such that for < e < Cq, 
the nonlinear system (4-3), (4-6) and (4-'^) has a unique solution 
satisfying the following estimate, 

< C({6-°-^|K}|hi + |{5^}|xo.r5xo..5(,) +6-l-^|{5^}U(,)}). 

(5.1) 

The proof of Theorem 5.1 is based on the existence of solutions to a 
linearized problem. Consider, 

< = + v\iy^ + A\^y + F\ e e fi^o, ^ > 0, 

(5.2) 




iALUe),r)-u^^'{-L^\e),T))^5\T), (5.3) 



u\C,0)=uliO, (5.4) 
with the compatibility condition 

Lemma 5.2. Assume that (5.2) satisfies PI and P2 in regular layers, 
and satisfies P3 and P4 in internal layers. Then there exists eo > 
such that for < e < eo, system (5.2)-(5.4) has a unique solution 
{u'}^^, u' e X^'\ni^ X R+,7), i e Z. Also 

I«ix2,i(-y) < C{e-o-5|K}|Hi +e-i-5|{F-'}U(^) + |{(5*}Uo.75xo.25(^)}. 

(5.5) 

The proof of Lemma 5.2 will be given at the end of this section. We 
now use Lemma 5.2 to prove Theorem 5.1. 
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Proof of Theorem 5.1. Let the solution in Lemma 5.2 be denoted by 

K} = ^(K}, {n). 

We use the Contraction Mapping Theorem to solve (4.3), (4.6) and 
(4.7). For {U'} e X2'i(7), let 

M=^(K}, m, {M\U\i,T)}). 

If |{f/^}|x2,i(^) < e''!, then from (4.5), \M'\x(^,) = o(e'i+i-5). From 
Lemma 5.2, if eo is small and < e < eo, |{'W*}|x2.i(7) < e^^- Thus 
maps a ball of radius e^^ into itself. 

From Lemma 5.2 and the second estimate of (4.5), it is also easy 
to show that JF is a contraction in such ball. Therefore there exists 
a unique fixes point for Using Lemma 5.2 again, we have 
(5.1). □ 

Proof of Lemma 5.2. There exist bounded extensions of Mq, to ^ G 
M. Without loss of generahty, assume that Uq G -ff^(M), -F* G X(M x 
M+,7). 

(I) First, consider (5.2) and (5.4) for ^ G M, and ignore the boundary 
conditions (5.3). The solution can be written as 

uHt) = e-^'X + r e^'^^-'^FUs)ds. 
Jo 

In regular layers, using Lemma 3.11 with ao — —0^/2, it is easy to see 

kix2.i(7) < CQuHh^ + |F^|x(7)). 
In an internal layer, we make the following spectral decompositions. 

ui = a'4>' + M^-^, 
F\s)=ff{s)ct>' + F^^{s), 
u\t)=u\{t)+uI{t). 

Here 0* is the eigenvector associated to the eigenvalue A*(e), u\{t) is in 
the eigenspace spanned by 4>\ the functions ^2(7"), u^^ and F*-'-(s) are 
in the spectral space corresponding to ReA < —a < 0. 



uUt) = e^' 



Jo 

Using Lemma 3.11 with ao — — a"/2, it is easy to see 

141x^.1(7) < C{\u},^\m + |i^'^k(7))- 

u\(r) = Q;^e^'(^)>^ + ( e^*^'^^^''^ pHs)ds)(l)' ^h + h 

Jo 
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Using the fact — - — < 7 < 0, we can verify that 



Let AP'{s) = (3'{s) - f3'{oo). Then e-^^A(3' G L\R+). 
Jo 



= ji + J2 + J3. 



A*(e) A*(e) ^ ^ ^ 



The first term is time independent with |Ji|ff2 < — |F*|x(-y). Since 

|e-7reV(e)r|^2 < I / 



< 



C 



l7l 



1.5 I 



k(7)- 



Since 



L2 



we have 



Therefore 



|wlU2.i(7) < C(e-°-Voli^i + 



Let the solution of (5.2) and (5.4) be denoted by u'' = u\ + u\. We 
have shown 

We now restrict to the domain = [— L^(e), L^(e)]. At the 
common boundary of and Ojj^ there is a jump 5* = ^ ('U*(L^(e))- 
u^+^(-Lj^^(e)). We can show that 

|5Vo.^«xo.25(^) < C(|4|^^i + |F^U(7))- (5.6) 
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In fact, in regular layers, using the Trace Theorem, 

Based on the Trace Theorem again, the above estimate is also valid for 
u\ in internal layers. For u\ in the internal layer, let 

u\(t) = h'(T)<lf. 

(4) = ■ 

Here |/i'(T)|iji(R+) < Ce~^-^{\uQ\H-^ -\-\F'^\x{'y)) ■ At the boundaries, |^| = 



1^^0.75 



+ |0||//o.2B < Ce-°^" ' < Cel•^ if e is sufficiently small. Thus 



<W\m{\<i>\±eP-')\ + W^{±e^-')\) 

J^0.7BX0.26 

< C(|<|^i + 

This proves (5.6), since — u\ + u\. 

(II) Consider a boundary value problem with zero initial condition 
and zero nonhomogeneous term, 

< = 4^ + ^'(04 + ^'(0< (5-7) 
u\Q) = 0, (5.8) 

5' = iu\LUe)) - u^^\-L^\e))) = 5^- S\ (5.9) 

It is obvious that 7ri5^(0) = 7ri5^(0) - 7ri5^(0) = 0. Also, due to (5.6), 

|5*|X0.75X0.25(^) < |5*|x0.75x0.25(^) + C ( | {Mq} | J^l + \{F'}\x{'y)) ■ 

(5.10) 

Let the solution to (5.7)-(5.9) be {u'}°°^. We want to show that 

|m*|x2,1(^) < C|{5'}|x0-'i'5x0.25(^). (5.11) 

Then from (5.10), and the fact -u* = + is a solution to (5.2)-(5.4), 
we will have proved the theorem. 

To prove (5.11), we write the solution as — u\ + ul. The function 
u\ satisfies a time independent system. 



n2n 
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The solution for that system uniquely exists and satisfies 

\u\\h^<C\{6\oo)}U2,.. (5.12) 

By converting the elliptic equation into a first order system of ODEs, 
the proof can be derived from [23], thus, will be skipped here. 
The function ul satisfies (5.7) with 



uiiC,0)^-u\{0. (5.14) 



I {<{LU^)) - ^n-L^'m = d\ (5.13) 



Here d'{T) = 6\t)-5\oo) E i/o.75xo.25(^^)_ ^5_;l4), the compatibility 
condition 

7r^d'{0) = -7Ti6'{oo) 

is clearly satisfied. 

To find ul, let ul = ul+ul- We want to satisfy the initial condition 
(5.14). Therefore, ul has zero initial condition, so that Lemma 4.1 

can be applied on {u\}. Assume that u\ has a bounded extension to 
H^(M.). Let Qq and Qg be the spectral projections associated to A*(e) 
and Re(A) < —a. Let 

ui = QquI + Q,ui 
Let Qoul{r) — hl{T)(j)'^ where 0' is the eigenvector, corresponding to 

A^(6). 

where |/i3(0)| < C|(5oWi|//i- Recall that the spectrum of A!' in "RQs is 
in Re A <—o. We have 

C 

|«3k2,i(7) < -^\Qqu\\h^ + CslQXIi^i- (5-15) 

Restrict u\ to the domain = [— L^(e), L^(e)]. At the boundaries, 
there is a small jump. 

We can show that 

141^0.75X0.25 < C|{wi}|jfl < C|{5X00)}|M2n. (5.16) 
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In fact, the traces of QsU\ at ±L^(e) is bounded by \QsU\\m. And 
using the exponential decay of 0* and 0| as |^| oo, the trace of Qoul 

|/i'3(T)0^(±e^-^)|j,o.75xo.25(^) < C|/i^|^i(|0^(±e^-i)| + |4(±e^-^)|) 

< c|(5oiti|i/i- 

Thus at ±I/J„ (e) , the traces of Qqu^ {^) are also bounded by C | {m\ } | //i . 
This proves (5.16). 

We now look for ul that satisfies (5.7), (5.8), with 

] {ul{r{e),T) - u^+\U^\e),T)) = d\T) - 4(r) = D\t). 
V^^/ (5.17) 

Observe that 

7ri(i^(0) = -TXi~5\oo), 

n,dm = uliLUe), 0) - ul^\-V^\e), 0) 
= -7ri5^(oo). 

Therefore, the compatibility condition 7ri£)*(0) = is vahd. Moreover, 

|-D*(t) I j:/0.75x0.25(-^-j < | (i* | j:/0.75 xO.25 (^-j + | (^3 | J|/0.75 xO.25 (--j,) 

Using Lemma 4.1 to (5.7), (5.8) and (5.17), we have a unique solution 
{u\} with 

< C|{-D*}|ho-75xo.25(^) (5.18) 

< |{(5'}|^0.76X0.25(-^). 

Recall that = -u*^ + ul, and = + M4. The part (II) of the proof 
is completed by combining (5.12), (5.15) and (5.18). □ 

6. Correction in Finite Intervals 

The main results of this section are Theorems 6.1 and 6.2. 
Recall from §5, in order to use Theorem 5.1 in [t/, 00), tf > log(l/(Coe^))/7. 
Define 

r = [log(l/(Coe^))/e7AT] + 1, 
At = eAr, tf = rAt, 

Ij = \jAt, {j + l)Ati < J < r - 1, 7oo = [tf, 00), (6-1) 
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where [x] is the integer part of x. As in §4, the coordinate change 
9^(0) ^'{1) maps X Ij to x The solutions of (1.1) in Ij 
arc equivalent to the solutions of (4.3), (4.6) and (4.7) in intervals I-l. 
More precisely, the solutions are indexed by j and will be denoted by 
u^j. The initial data, jumps and forcing terms are Mj(0), S^, and 
the coefficients are ^j(0> ^*(0- The domain for equation (4.3) is 
n] X P^. Recall that = [0, At], < j < r - 1, and /~ = 1R+, and 
Q] = (-L;.(e),L;.(e)) where L;.(e) = Af{t)/2e, t^jAt. 
In the new notations, consider, 

4r = + + ^KO^ + -A/JH, e), e e ^, r e 

(6.2) 

[d^ ^' ^) - ^^?'(-i^r(e), r, 6)] = 5i(r), (6.3) 

w}+i(e, 0) = t.;.(6, At), ..^(e, 0) = ul^iO- (6.4) 

In (6.4), it is understood that -u* = u,^, and (.C, r) {^i,Ti) is the 
change of coordinates 1K*(1) ^ 9i(0) W(l) that maps 



With t = jAt, define. 



z/1 : ,T = m*(t + At) +e^, 
z/2 : y = X — S(a;, t + At, t, i), 

^^3: ei = (z/-rn'(^))A- 



We have 



tx}+i(e,0) = M}(z/3oz/2oz/i(0,Ar). 

We have showed, in §4, that (4.3) satisfies PI and P2 in regular 
layers, and satisfies P3 and P4 in internal layers. Also satisfies 
(4.4). In each finite time interval, the following existence of solutions 
holds. 

Theorem 6.1. For any At > 0, let the intervals Ij, < j < r — 1, 

be constructed by (6.1). Consider (6.2) and (6.3) in an time interval 
— the image of Ij by the change of variable r = {t—jAt)/e. Assume 
that g] e L'^ip^j x l(), u){Q) G H\n)), and 5] e //o.75xo.25(ji)^ ^^^f^ 

\{g]}y = o(e), \{ui{Q)}\m = o(e), 1^0.75.25 = o(e). Then there 
exists eo > such that for Q < e < cq, (6.2)-(6.3) has a unique solution 

{u]}^^, u) G X li). Also 

\u)\H-^ < Ce'^^^AT{\{u'j{0)}\m + li^jllL^) + C|{5i}|^o.75xo.25(4). 
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In the time scale t, the length of the interval [0, At] is only eAr — 
0(e). Since ?i*(Ar), after a near identity change of variable, is the 
initial data for the next time interval, and since r — oo as e — > 0, the 
result in Theorem 6.1 is not accurate enough to guarantee the existence 
of solutions in I^x,. We will prove a sharper estimate in Lemma 6.4 
showing that 

m(AT)|Hi < (1 + Ce)|{4(0)}Ui + . . . . (6.5) 

With this, even if r = 0(- log(e)/e), using (1 + Ce)^/^" < e. the 
accumulation of error will not grow too fast as e — > 0. We show in 
Theorem 6.2, that the accumulation error at = rAt is small enough 
so that the existence of solutions in Ij is guaranteed by Theorem 6.1, 
and in is guaranteed by Theorem 5.1. To this end, let us recon- 
sider (6.2)-(6.4) with an emphasis on the norm |ti*(AT)|iji in term of 

|{4(o)}|ffi. 

Theorem 6.2. There exist At > such that if Ij and loo o-f^ defined 
by (6.1) using this At, and then we have the following properties con- 
cerning the solutions of (6.2)-(6.4)- If ri > 1.5, Ji = where C is 
from (6.5) and from Theorem 5.1, |{mo(0)}|//i = o{e^^^^-^^'^^), |{5*}|j:^o.75xo.2b-|- 
\{g]}\L2 = o(e''i+i-5+^i), then system (6.2)-(6.4), < j <r - 1, has a 
unique solution for < e < Cq, where eg > is a small constant. The 
solution in each time interval satisfies 

In particular, |{Mj^(0)}|iji = o(e''^+°-^). 

The proofs of Theorems 6.1 and 6.2 are based on Lemmas 6.3 and 
6.4. We present these lemmas first. 

Since all < J < r — 1, are identical, they are denoted by for 
simplicity. Assume At > 1, in = [0, At]. Consider the following 
initial value problem, 

= u^^ + V{i)u^ + A{i)u + F(e, t), e e t e 7„ 

(6.6) 

^(e,0) = «o(e)- (6.7) 

Lemma 6.3. Assume that the linear operator A : IP' L'^, u — > 
u^^ + V{^)u^ + A{^)u has a simple eigenvalue Aq that is close to zero, 
all the other spectra are contained in {Re\ < —a} for some constant 
a- > 0. Let -a < ^ < 0, and n = max(0, iieAo). Let F e L^(R x 7^) 
and uq e 77^ (R). Then the solution of (6.6), (6.7) satisfies 

k|i/2.i(Kx/o < Ce^"^' AT{\u^\m + IT^Il^). (6.8) 
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Moreover, let Qo and Qg he the spectral projection corresponding to Aq 

and the spectra in {Re\ < — o"}. Since u : Ir L^(R) is continuous, 
we can write u(t) = h{r)(f) + Qsu{t), where h{T)(p — Qqu{t). Then 

(a) \h{T)\ < e^'^^^\h{0)\ + e^'^^/¥\F\L2). 

(h) \Qsu{t)\h^ < C{e^'^\QMm + W^)- 

Proof. Let uq = h{0)(t) + QsUq and F{() = p{()(f) + ^^(C) be the spec- 
tral decompositions corresponding to Qo and Qg. Applying Qs to the 
variation of constant formula, we have 

Jo 

When restricted to TZQs, ^ is a sectorial operator with Rea{A) < 
—a < 0. Thus, e-^'^ is exponentially stable on TZQs- Therefore, we have 
(b). And also 

\Qsu\H2,i(Rxir) < C{\uo\h^ + \F\l2). (6.9) 
Using Qo to the variation of constant formula, we have 

h(r) = e^°^h(0) + r e^°^^-^^ mQdC- 
Jo 

Jo Jo Jo 

< e^^^V^\(3\Lm^) 

< Ce''^V^\F\L.(j^xi^y 

Observe also that \e^°^h{0)\ < e^^^°^|/i(0)|. We have proved (a). 
Based on (a), it is elementary to show 

\h\HHu) < Ce^^^(x/A^|/i(0)| + At\F\l2). (6.10) 

Since \h{0)\ < C\uo\h^ and iPli^ir) < C'I-^|l2(mx/,), 

\h4>\mA(Rxir) < C\h\Hi{i,)\4>\m 

< Ce''^'-{VA^\uo\H^ + At\F\l2). 

Recall that At > 1. Thus 

m + |-^|l2)- 

(6.8) follows from the above and (6.9). □ 
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Consider a linear initial boundary value problem in one of the finite 
time intervals. Let Ir = [0, At] and = (-Li(e), Lj(e)). 

(6.11) 

[d^ - u^^\-V^\e))) = S\ (6.12) 

^U0) = ii^(O, (6.13) 
with the compatibility condition, 

Assume that the coefficients of (6.11) have been extended to ^ e M 
by constants. Let the operator on -u* defined by the right hand side 
of (6.11) be denoted by A' : L'^{R) L'^{R). Let e H\R), F' e 
L\R X Q and 5' e //0-^^^o-25(/^). 

Assume that (6.11) satisfies PI and P2 in regular layers and satisfies 
P3 and P4 in internal layers. In internal layers, for any u e H^{M.), 
let u = QqU + QsU be the spectral decomposition. Let QqU = hcjf 
where 0* is the eigenvector corresponding to A*(e). Then \h\ + \Qsu\h^ 
is an equivalent norm io \u\hi. In regular layers, for convenience, let 
QsU ~ u, h — 0. Formally, we still have the above equivalent norm. 
Let /I — max{0, supj[ReA*(e)]}. 

Lemma 6.4. There exists eo > such that for < e < cq, (6.11)- 
(6.13) has a unique solution {u^}°^ooi ^ -ff^'^(^- x 1^) for all i e Z. 
Also 

(6.14) 

Moreover, i/ At is sufficiently large, then the solution admits a bounded 
extension to ^ e M such that if QqU^ — then, 

|/i^(AT)| + |gy(AT)|j,i(M) < (l + Ce)sup(|/i^(0)| + |gy(0)|^i) 

' (6.15) 

+ C\{5'}\ ^0.76X0.25(7^) + C\{F'}\ 
The constant C depends on At, hut not e. 

Proof. The solution will be written as — -\-u^. Using the spectral 
projections, 

u\T) = h\T)cl>' + QsU\T), 
U\T)^h\T)(t)' + QsU\T). 
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Let u^{t) — e-^^'^Mo + e'^^'^'^~^^ F^{s)ds. In internal layers, using 
Lemma 6.3, (6.8), we have 

|wVm(mx/o < Ce^''^AT{\ul\Hi + \F%2). (6.16) 

(6.16) is also valid in regular layers. 

In internal layers, u^{t) — ff{T)(f)^ + QsU^{t). Note that ^'(O) = 
Uq, ^*(0) = /i*(0), QsU^{0) — QsUq. We have, using Lemma 6.3 (a) and 
(b), for some —a /A < 7 < 0, 

\h'{AT)\ < e>^^^mO)\ + VA^\F%.), 

\Q,u\At)\h^ < C[e^^^|QXI/fi + \F%A- 

(6.17) is also valid in regular layers, with — 0. 
At^ = ±L}(e), 

\u\±me),r)\ < \¥{rWi±me))\ + \Qsu\±me),r)\. 
It is clear that from the Trace Theorem, 

I i^D^ Qs'U'(±iv}(e), •)|i/0.7Bx0.26(^) < C\QsU^\h2A{^) 

<C(|QX|^i(M)+el^l^^|F'U2). 

The second inequality is based on Lemma 3.11. 

Though Qqu'^ does not decay exponentially in time, we can still con- 
sider weighted norm in the finite interval 1-^. Observe that 

Using (6.10) and |0*(C)l + l^€0*(C)l ^ Ce'^^l^l for some 71 < 0, we have 

I Qou\±me),-)\HO-rsxo.2^^) < Ce(^+l^l)^^e^^^'"'(VA7|/iX0)|+Ar|F*U2). 

Let eo > be small so that Ce^^'+^^^^^^+^''^~' At < e, whenever < 
e < eo. Then the above is bounded by e(|/i'(0)| + |F*|i2). Whence, 

I MX±^}(e)'-)Uo.T5xo.25(^) < e\h\0)\+C\QsUl\m+Ce\''\^^\F%2. 

The above estimate is also valid in regular layers, with /i*(0) = 0. 
Define, 



^ - [d^) {u\L){e))-u^^\-Lf\e))). 
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Consider a new boundary value problem with zero initial condition 

= Au , 




= 0. 



Since It is a finite interval, 5* is also a point in H^-'^^^^-'^^{It-, 7), with 

|(5*|^0.76X0.25(/^^^) < Ce'''^''^^|(5*|//0. 75X0.26(7^). 
\5 |f70.76x0.25(/^^^) < Ce'''''^^|5'|ffO. 75X0.26(7^) 

+ sup{e\h\0)\+C\Qsui\H^+Ce\'^\''^\F%2}. 

i 

We now use Lemma 4.1, case(2), with 7 < 0. If e is small, the system 
has a unique solution {w^j^fj^, with 

I^^1if2,l(nix7r,7) - }|i70.76x0.25(^). 

First let 7 = in the above. Combining that with (6.16), and 
recaUing — + 11,^, we have (6.14). 

Now let — (t/4 < 7 < again. Using a bounded extension of T to 
^ e M. At r = At, we have 

|^'(At)| + \Qsu\At)\hhr) < C\u\At)\hi < Ce'^^-|S%2,i(nix7.,7) 

< C(|{5^}|770.75xo.25(7,) + e'^^^sup{e|/i^(0)| + \Qs4\m + e'^l^^jF 

i 

Combining this with (6.17), and recalling that — u'^ + u^, we have 
\h'{Ar)\ + \Qsu'{Ar)\m < C|{5^}|770.76xo.25(7,) 

+ sup{Cie^^^\Qsui\Hi + {l + Ce)\h\0)\ + Ce^^^VA^\F%2}. 

i 

Here, we used the fact |;Li| < Ce and e^'^'^ < 1 + e^^^fiAr to simplify 
(6.17). Observe that Ci does not depend on e or At. Let At be 
sufficiently large so that Cie^^'^ < 1. Prom this, estimate (6.15) follows. 

□ 

Proof of Theorem 6.1. Let the solution of Lemma 6.4 be denoted 

by 

Let W e H^'^{fl' X It) with \W\h2.i < e. Let 

W =.F(K}, m, {Af\U\U+er,e)}). 
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Based on (4.4), |jV*|l2 = o(e). Thus, if e is small, from Lemma 6.4, 

|^'*|i/2,i < e. Therefore, jFmaps an e-ball in H"^'^ into itself. From (4.4), 
it is also clear that JF is a contraction, provided that e is small. Thus 
there exits a unique fixed point {u''}'^^ for J^. The desired estimate 
follows from Lemma 6.4 and the fact 

< \9%2+0{e)\u%2^ 

□ 

Proof of Theorem 6.2. Let At be as in Lemma 6.4, and let Ij, < 
j < r — 1, /oo be defined by (6.1) from this Ar. 

The spectral projections used in Lemma 6.4 actually depend on the 
layer index i and the index j of time intervals, and shall be denoted 
by Qo"' and Ql'^ ■ In regular layers, of course, Qq-' = 0. Let the critical 
eigenvalue be A*(e). The eigenvector 0* corresponding to A*(e) and the 
covector ipj associated to the adjoint equation are normalized so that 
|0j|L2 = 1 and < >— 1. The perturbation theory of eigenvalues 

and eigenvectors yields, in internal layers, 

10} - + iVi - V'i+iU^ = 0(e). (6.18) 

Consider the solution of (6.2)-(6.3) as in Theorem 6.1. When |{w*(0)}|^fi = 
o(e), |{5*}|j^o.75xo.25(7^) = o(e) and |{5']}|l2 = o{e), the unique solution 
in that theorem also satisfies |{ii*}|j^2,i = o(e). Now using Lemma 6.4, 
with f/j replacing F*, we have 

|/i}(AT)| + |g:'%}(AT)|^i(M) < (l + Ce)sup(|/i5(0)| + |g:'X(0)|Hi) 

i 

+ C|{(5}}|^^0.r5x0.25(,^)+C|{A/J}U2. 

But from (4.4), 

< \{9]}\l2 + Ce{\{u]m\H^ + \{g;}\,. + |{<5]}|ho.75xo.25(,^)). 
Therefore 

|/i}(AT)| + |g:'X.(AT)|Hi(R) < (n-C6)sup(|/i}(o)| + |g:'X(o)|//i) 

i 

+ C'|{5'}}|l2 + C|{5]}|ij0.75x0.26(7^). 

Consider now the initial data in the next time interval [{j+l)At, {j+ 
2) At], 

u]^,{^,0) = ^-(6, Ar) = u]{us o U2 o At). 

Cf. (6.4). From the definitions of the mappings ui, 1^2, i^3, we can verify 
that 
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where C is a smooth function satisfying 

\aeO\<Ce; \Cy<C. (6.19) 

The constants C depend on At but not on e. The proof of (6.19) uses 
the coordinate change formula in §4 and is tedious,. Details will not 
be given here. Based on (6.19), we can show, 

|^(6)-^(0l<Ce|e||%|oo. (6.20) 

For any g e H^{M.), since ^ — > .^i is a near identity mapping, we can 
verify that, 

|gr o 1/3 O 1/2 O Ui\hI < (1 + C6)\g\Hi- (6.21) 
In regular layers, from (6.21), we have, since Qo''^^ = 0) 

In internal layers, observe that for g e L^, Qo^~^^g — (/f^ '0j+i5'c^O^'^}+i- 
Thus, 

hi^M = |v'i+i(0[^K^^)4(^i) + ^^S(^l'^^)]^^■ 
/V'i+l(O</'5(^l)^^=/V'i(O0}(Oc^e 

+ /(^;+i(c) - mMiOd^ + /^i+i(0}(ei) - <p]{o)d^ 

< 1 + Ce. 

In proving the above, (6.18), (6.20) are used. 

Observing that the up to the second derivatives of 0}(C) approach 
exponentially as ^ — > ±00, uniformly with respect to i, j. Prom (6.19), 
and the Mean Value Theorem, we have 

mi) - m)\ + \9dmi) - <i>m\ ^ mie-^^^^^, 

for some 71, X > 0. Therefore, 

10* O Z/3 O Z/2 O Z/i - (f)'j\Hl < Ce. 

Using the above, evaluating Q'^g''^^^] — Qs''~^^{4>) ~ ^j+i)' then using 
(6.18), we can prove that 

IQ:'''''(0}°^3oz/2oz/i)|hi <Ce. 
One can similarly show that, 

I / vi+i(OQ:'''^}(ei,AT)dei <^^e|g:'^\.5(-,AT)ii.. 

m'''\Q'j'u}{us o o Ar))|^i < (1 + Ce)|g:'%}(-, At)|^i. 
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Combining all these estimates, we have shown, 
( l^;+i(0)l WA + Ce Ce \( \h]{AT)\ \ 

Therefore, 

l^}+i(0)| + |g:'^-+^4+i(0)|^i < (1 + Ce)i\hiiAr)\ + |Q>}(Ar)|HO- 
Combining the above with Lemma 6.4, we have 

< (i+Ce)sup{|/i5(o)| + |g>}(o)|Hi}+C7(|{5i}|^o.75xo..5 + |{(?;.}UO- 

i 

By induction, the initial condition {u^ (0)} for the j'-th time interval 

satisfies 

\h]m + |g:'%}(o)Ui < (1 + cey sup{\him + iq^u^^m 

i 

k=0 

<{l + Ceysup{\hi{0)\ + \Qfui\m} 

i 

+ ^^^^(imi/^o.™ + iMi^.). 

Since j<r< log(l/(Coe2))/(e7Ar) + 1, 

(1 + Cey < (1 + Ce)(l + Ce)^e^''^^^/iCoe'))^ 

< (1 + C6)e^°^(^/^^-^»* 

< (1 + Ce)[Coe^]^ < Ce'^. 

The estimate in Theorem 6.2 follows easily. □ 

7. Application to Singularly Perturbed Reaction 
Diffusion Systems 

We first review the construction of matched asymptotic expansions 
for a wave- front-like solution for system (1.1), cf. [25]. We then show 
that the Spatial Shadowing Lemma can be used to obtain exact solu- 
tions near the formal expansion. 

As in §1, assume f{u,x,e) = YfjL^^^ fjiu^x), and f{u^x + Xp,e) — 
f{u, X, e) for some Xp > 0. We will consider internal and regular layers 
occurring alternatively in the x direction. 

Assume there is a partition R = [J'^_^[x^, x^~^^] that is periodic with 
respect to i, compatible with the period of /. That is, there exists ip 
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such that — X(,-\- Xp. A C°° function p'^{x) is defined on x^\ 
with /o(p*(x), x) = 0. Also assume that + Xp) — p^{x), x e 

[x^~^, x^]. 

H* 1 Rea{fou{pKx),x)} < for x e [x'-\x'], i e Z. 
Using H*l, a formal expansion 

oo 


for (1.1) can easily be obtained. Note that the expansion is time- 
independent, since f is so. 

Using a stretched variable ^ = , we assume that the 0th order 

e 

expansion of Eq. (1.1), 

u^^ + fo{u,x')^0. (7.1) 

has a heterochnic solution g*(^) connecting p*(a;') to p^^^{x''). Assume 
H* 2 The hnear homogeneous equation 

0ec + /o«(g*(O,^O0 = o, 

has a unique bounded solution q^iC,), up to constant multiples. 
From H*2, we infer that the adjoint equation 

V'«+/o"„(gi(o,^> = o 

has a unique bounded solution '?/'j(0 constant multiples. More- 

over, ^ exponentially as |^| oo. 

The following assumption implies that the heterochnic solution breaks 
as X moves away from x\ 

Under the hypotheses H*l to H*3, we can construct formal expan- 
sions for the stationary wave front positions and solutions in internal 
layers 

oo 

t Xj, Xq — X , 



oo 



The formal solutions in the internal and regular layers match at their 
common boundaries. See [23, 25]. 

When x is in a neighborhood of x\ we look for a travelling wave 
solution with wave speed V'^{x). It is also of interest to find out 
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conditions to ensure the wave front to move towards x^. Consider 
ji : l2(^) ^ ^2(^)^ defined as 

H* 4 A = is a simple eigenvalue for i e Z. There exists a > 
such that 

a{A') C {A = 0} U {ReA < -a}. 

Hypothesis H*4 implies that q^C) is a stable equilibrium for 

Mr = + /o(M,a;'), (7.2) 

modulo a phase shift in ^, see [9, 10, 11, 12]. Since A = is simple, we 
have 

/oo 
mmodi + 0. (7.3) 
-00 

H* 5 J(xO = [roo^r(04(0^e]-'/-°°oo^^^(0/o.(g^(0,^O^e > o, ^ g z. 

Condition H*5 implies that the wave front is moving towards a;*, if 
X is near x^ . In fact, it is shown in [25], 

^ = -/(.■). 

ox 

An important consequence of (7.3) is that if x is near x*, there exists 
a unique V = V'^{x) such that (7.2) has a travelling wave solution 
q^{C,x), with the wave speed V, satisfying 

u^^ + V'{x)u^ + fo{u,x)^0. (7.4) 

The function g'(^,x) connects p^{x) and p^'^^{x) as ^ — > ±00. More 
precisely, under the Hypotheses H*l to H*5, there exists an open in- 
terval O*, containing such that the following properties holds. For a 
proof, sec [25], Lemma 4.1. 

(i) p*(x) and p'^^{x) can both be extended smoothly to with 

fo{p'{x),x) = 0, 

^Ga{fou{p'{x),x)) < 0, j = i,i + l. 

(ii) There exists a C°° function : O* — * R, such that for each 
X G O*, Eq. (7.2) has a unique heteroclinic solution q^{^, x), connecting 
p^{x) to p'"*"^ (x) , with q^{Q,x) — g*(0) ± 5^(0). In particular, V'^{x^) = 
and q^{^,x^) — q^{^). Moreover q^{^,x) is C°° in both variables. 

(iii) The linear equation 

<f>^^ + V'{x)<l>^ + fou{q% x'),x)<f> = 0, 
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has a unique bounded solution up to constant multiples. The 

adjoint equation 

- v'{x)i;^ + f^^{q% x'),x)^ = 0, 
has a unique bounded solution ipi{^,x), normalized so that 

/oo 
-oo 

x) is a C°° function of x in the space C(]R). 

(iv) The densely defined closed operator A^{x) : — > L^, 

A\x)u = u^^ + V\x)u^ + fou{q\^, x),x)u, 

has A = as a simple eigenvalue with eigenvector x). All the rest 
of the spectrum are contained in ReA < — ckq for some ckq > 0. 

(v) For xeO\ieZ, 

/oo 
''PI{^,x)fouiq\^,x),x)d^>0, 
-oo 

As a consequence of (v), V^{x) > 0, (< 0) if x < x\ (> x'), and 

^^^J^ = -l{x), xeO\zeZ. 
ox 

Under the Hypotheses H*l to H*5, we can construct formal series 
solutions and wave positions for (1.1) in the i-th internal layer, 

oo oo 


The formal expansions are to be truncated to the mth order, as in 
(1.2), where Uj^{x,t) = u^{x) is, in fact, t-independent. Based on 
these expansions, we define a formal approximation {w*} by (1.5). In 
the rest of the section, we show that Hypotheses H1-H4 in §2 are 
satisfied by this {tf*}. Besides, by choosing large m, the results in 
[25] indicate that j2, as in (2.4), can be arbitrarily large. Therefore, 
the existence of exact solutions near {w*} is ensured by the Spatial 
Shadowing Lemma as in §2. 

To the zero-th order, the wave position obeys the equation 

d^ Mil \ 

Tt = 

with an initial condition a;(0) = ?7o(0) ^ x — x'^ \s di, stable equihb- 
rium due to the condition Dj.^ (x*) < 0. Therefore 

mt) - x'\ < Ce-^\ 
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for some C, 7 > 0. In fact, it is shown in [25] that 

\Dtrt,{t)\ + ^{t) - r/;(oo)| < Qe"^*, j G Z, 
for some Cj > 0. Moreover, in the weighted norm, 

sup{(l + I uf{i, t) - uf{i, 00)} < C,e-^\ 

for some j > 0. In the domain |^| < e^~^, we have (e|^|)-^ -C 1, and 

m 

^e>f (e,T) -lif (e,oo)| < Crnc-^K (7.5) 

j=Q 

Also 

m m 

J2e^D,r^}{t)\+J2e^\v]{t)-v]M\ < Cme-"^'. (7.6) 

j=0 j=0 

The constant Cm depends on the order of truncation m. It folfows from 
(7.5) and (7.6) that HI in §1 is satisfied in internal layers. HI is cer- 
tainly satisfied in regular layers where the expansions are t- independent. 

Prom H*l, the Hypothesis H2 in §2 is satisfied if ty* = p''-{x). In 
general, adding higher order terms introduces a perturbation of 0(e). 
From the perturbation theory of eigenvalues of the matrix /ou, H2 in 
§2 is valid if is an m-th order truncation of the formal solution. 

Let XltT ^''^f^ii-i t) be the inner expansion of the outer solution t^uf'^{: 
[25, 7, 8]. For any > 0, if e is small, we have, in the region |,^| < e^~^, 

rn rn 



= 0(e^("+^)) < ^/3. 

It is shown in [25], based on the matching of expansions in internal 
and regular layers, there exists 71 > such that 



-71 1^1 



In the domain |^| < e^~^, e^{l + |^p) <^ 1, thus the above is bounded 
by Ce~'''il^l. For > 0, there exists > such that the above is 
bounded by for < ^ < e^~^. Observe that when |^| < e^~^, and 
e is small, 

m m 

I Y: e^ufiv'it, e) + t) - ^ e^uf{v\t, e) + e^, t) \ 


< C|e^ - e^l < At/3. 
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Combining the three estimates, we have proved H3, §2 for €^~^ > ^ > 
N. The proof for — e^~^ < ^ < —N is similar. 

It remains to show H4 of §2. From the Property (iv) of 0\ when 
e = 0, A^{t) has a simple eigenvalue A = 0, all the other spectra 
are contained in ReA < — ctq. By a standard perturbation analysis, 
for any 5 > 0, there exists eo such that for < e < eo, A^{t) has a 
simple eigenvalue |A*(e)| < S, while all the other spectra are contained 
in ReA < —ao + S. Moreover, the eigenvector corresponding to the 
critical eigenvalue can be written as + z, with z{0) _L g^(0) and 
\z\ < 5. it is also clear that A*(t, e) = eXl{t) + O(e^) due to the smooth 
dependence of eigenvahic on e. Details will not be rendered here. We 
now show the following important result. 

Lemma 7.1. 

Proof. The eigenvalue A and eigenvector 2; + g| satisfy 

X{z + ql) = {z + ql)^^ + Dtri\t, e)iz + + Uu\ + 6^, e)iz + g^). 

Here — e^Uj and 77* = e^rjj are formal expansions for the 
solutions and wave front in the ith internal layer. Using the fact is 
an eigenvector corresponding to A = 0, we rewrite the equation as a 
nonhomogeneous equation for z with forcing term that will be zero if 
e = A = 0. The operator 

is Fredholm, in the space L^(l^l) = {u : u{^)/{l + \^\) G L^} of weighted 
functions, with one dimensional kernel and co-kernel. The reason 
to use weighted norm is that the partial derivative of the right hand 
side with respect to e yields a factor ^. Lyapunov-Schmidt method can 
be used to deduce that there exists a unique solution z{{)) ± q\{^) if 
a bifurcation equation G{\ e) = is satisfied. The partial derivatives 
can be obtained as Melnikov type integrals. 



^(^(0,0) 

d\ 
^^(0,0) 

de 



/oo 
-00 
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Introducing a variable ^, we can rewrite the above as, {x = rjoit)), 
+ fo.(q\^ + ^, x),x)(4(t) + ^ + + fi(q'(^ + ^ x),x)}d^ 

/oo _ _ _ _ 

^li^ + [fou{q'{^ + ^, x),x)u\^{^ + + fox{q\^ + ^, x),xM 
-oo 

Observe that Uq = q\ and 
<t = + Dtvi{t)u\^ + fouu\ + Dtrj{{t)ui^ + /o. • {4{t) + + fu 

which can be obtained from expanding (1.1) in powers of e, see [25]. 
We have, 

dG 



de 



/oo 
-oo 

Observe that H = + DtrfQ{t){u\^)^ + fouU\^ is in the range of 

the operator (when e = 0). Thus J^^tjjlHd^ — 0. Also observe 
— Dtq^{^,rf{t)) — ql.V'^{x). Therefore, we have 

The desired result follows from the following formula, see [25], 



/OO fOO POO 

i^IqidC]-'{[ / i^JqUdC]v\x) - / vr/o.rfe} 
-oo J—oo J— CO 



/oo poo 
-oo J —oo 



dx 

□ 

Corollary 7.2. // ^^g^^ ^ 7^ 0, then there exists a neighborhood of x^ 
where 

dVUx"-) 
sign{\l{t)) = sign{ — — — ) 

The proof of Corollary 7.2 is based on = 0, and V^{x) is small 

if X is near x\ From H*5, it is clear that Aq(oo) < 0, since r]Q{oo) = x\ 
The method of introducing a new variable ^ in the proof of Lemma 7.1 
has been used in [23]. 
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